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stood to apply SAR image spectra with confidence for wave studies if supported 
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A theoretical framework is developed for the assimilation of arbitrary wave 
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Abstract 


A study on the applicability of ERS-1 wind and wave data for wave models is carried out 
using the WAM third generation wave model and SEASAT altimeter, scatterometer and SAR 
data. A series of global wave hindcastsare made using as input driving fields. 


(i) the six-hourly averaged surface stress fields derived by Atlas et al (1 987) for the full 96- 
day SEASAT period 7 July - 1 0 October 1978 by assimilation of scatterometer data with 
conventional meteorological data using the Goddard Laboratory of Atmospheres 4 : x 5 
9 level atmospheric model, 

(ii) the 1 000 mb wind field from the same GLA analysis for August, 1978, 

(iii) the 1000 mb wir >d field derived by Anderson et al (1987) for the period September 6-17, 
1978, using more sophisticated assimilation techniques and the higher resolution T63 

(1 875’ X 1 187°, 19 layer) model of the European Centre for Medium Range Weather 
Forecasts, and 

viv) a subjectively analysed scatterometer wind field produced by Woiceshyn et al (1987) for 
the period September 6 • 20 

The four hmdcasts are intercompared and verified against altimeter wave height and 
wave buoy data In the northern hemisphere the hmdcasts agree reasonably well with each other 
and with the observations However, strong deviations between the hmdcast for the GLA stress 
field and the other three hmdcasts are found in the high wind belts in the southern hemisphere 
This is tentatively attributed to the six-hourly averaging applied in the computation of the 
surface stress fields 

It is concluded that wave models provide a very sensitive mechanism for identifying 
inconsistencies in wind field analyses The simultaneous operation of a wave model in 
conjunction with an atmospheric mode' in a joint wind and wave data assimilation system 
therefore provides a valuable data validation tool 

Comparisons are also made between SEASAT sar image spectra and theoretical SAR spectra 
derived from the hindcast wavespectra by Monte Carlo simulations Good overall agreement is 
found for 32 cases representing a wide variety of wave conditions in most cases the input wave 
spectrum is strongly distorted in the SAR spectrum by nonlinear motion effects T he linear 
treatment of motion effects (in terms of a velocity bunching transfer function) was found to yield 
unrealistic monotomcally increasing wave spectra in all cases (including swell fields) It is 
concluded that SAR wave imaging is sufficiently well understood to apply SAR image spectra 


with confidence for wave studies However, SAR data can be meaningfully interpreted only in 
conjunction with a realistic wave model and detailed computations of the mapping of the two- 
dimensional ocean wave spectrum into the SAR image spectrum. 

A new, closed nonlinear integral expression for this spectral mapping relation is derived 
which avoids the inherent statistical errors of the Monte Carlo approach The method may aiso 
prove to be more efficient numerically This is an important consideration in view of the 
extensive computations required for the analysis of the global SAR spectral data set produced by 
the ERS-1 AMI instrument operating in the wind/wave mode 

A theoretical framework is developed for the simultaneous assimilation of arbitrary wave 
data le g. wave heights, buoy spectra or SAR image spectra) in numerical weather prediction ana 
wave models Both wind and wave fields are modified simultaneously m accordance with the 
constraints imposed by the wave model The explicit integration of the adjoint wave equation 
required in the general formulation of the problem is avoided by using an approximate 
(diagonal) Green function, the elements of which are already computed as part of the wave 
model's implicit integration scheme 

An example of a simplified wave data assimilation scheme is presented in which only the 
wave field is modified The assimilated wave data are the global SEASAT altimeter wave heights 
for August, 1978 A considerable improvement in the wave forecast is achieved m the tropics and 
northern hemisphere, where most of the wave energy consists of swell However, the errors are 
only partially reduced in the principal generating regions of the higher latitude southern 
hemisphere, where the unmodified surface stress field regenerates incorrect windseas These 
deficiencies would presumably be remedied in a joint wind/wave assimilation scheme 


Introduction 


1 1 Background 

The present study contract represents a complementary investigation to a 
previous ESA two-part study on the feasibility of developing and implementing a 
comprenensive data assimilation system for the near real time analysis of ERS-1 wind 
and wave data It should therefore be viewed in the context of these study contracts 
(ESRlN Contract No 6297/86/H G E - 1 (SC) , Part I, Anderson et al., 1987, Part II, Janssen et 
a!., 1988) The work carried out in these three studies must also be seen <n a broader 
perspective within the framework of a longer term international programme to use 
EPS- 1 wind and wave data, together with conventional weather network data and 
other m situ data, to compute the global fields of all fluxes at the air/sea interface on 
a continual, operational basis, beginning 1991 These activities are coordinated m the 
Global Assimilation Programme for Air/Sea Fluxes by the JSC/CCCO Working Group 
on Air/Sea c luxes (cf Appendix) In the following a brief summary is given of earlier 
work which has provided the basis for the studies earned cut m the three ESA 
contracts and for the planning of the Global Assimilation Programme for Air, Sea 
Fiuxes 

1.2 Need for a data assim i lation system 

The extensive oceanographic surface data which will be provided by the suite of 
sensors on ERS- 1 open exciting new opportunities to tne oceanographic community 
But the preparation for the effective use of these data also poses a severe challenge 
This is vs/el! illustrated by SEASAT The wave hmdeasts presented m this study and in 
the companion study contract (Part I!, P Janssen et ai , 1 988) represent the first 
application of SEASAT data to a global ocean wave hmdeast — fully nine years after 
the launch of SEASAT! The computation of the surface stress fields used in these 
hindcasts for the 96-day period during which the SEASAT scatterometer was 
operating was only completed, with considerable effort and perseverance, eight 
years after launch (Atlas et al , 1987) The analysis was based on an off-line data 
assimilation system using the non-operational, relatively low resolution GIA 
(Goddard Laboratory of Atmospheres) atmospheric model r hese aata processing 
problems would have been still more apparent if SEASAT had flown, asoianned. for 
several years rather than three months Although SEASAT provided an impressive 
demonstration of the great scientific and operational potential of ocean satellites, >t 
also clearly revealed the problems the scientific community faces in using satellite 
data if provisions are not made for the timely implements ion of an effective 


operational data assimilation system capable of processing all incoming satellite and 
conventional data in a single quasi real-time operation 


It should be stressed that the mam delays in the production of the level III final 
SEASAT data set were not caused by the processing of the raw data or the provision 
of appropriate algorithms, but by the last, most complex element of the data 
processing chain: the data assimilation system The assimilation systems used by Atlas 
et al (1 987) and Anderson et al . ( 1 987) can furthermore still be viewed only as a first 
step ideally, such a system should combine all the wind and wave data provided by 
the satellite with all relevant data from other sources to construct simultaneously 
global fields of the surface wind and sea state. Both m the GLA system and in the 
more sophisticated ECMWF (European Centre for Medium Range Weather Forecasts) 
data assimilation system, wave information was not used in constructing the wind 
held Another important feature of an effective assimilation system is that it should 
be implemented in an operational global wind and wave forecasting environment. 

♦Although evidence of deficiences in the SEASAT scatterometer data had already been 
pointed out earlier (Woiceshyn et al , 1987), the present study, together with the 
previous ESA assimilation studies, revealed many further sources of error and data 
inconsistencies These can be routinely identified and corrected only by carrying out a 
complete data assimilation and data validation cycle within a quasi -operational 
setting, using state-of-the-art, high resolution atmospheric circulation and global 
wave models. 

1 3 Previous work 


The scientific community has long been aware of the potential but also the 
difficulties of implementing an effective ocean satellite data assimilation system 
Accordingly, since the early conception of ERS- 1 it has been preparing a programme 
and some of the necessary tools for the development and implementation of a 
comprehensive, quasi real-time, operational data assimilation system for the satellite 

The first step in this programme was the development of a realistic global wave 
model While the ERS- 1 community had access to the world's foremost global atmo- 
spheric model and atmospheric data assimilation system at the ECMWF, a global 
wave model of similar capability did not exist. The only global wave models in 
operation in the early eighties were first or second generation models, whose basic 
shortcomings have been identified in the Sea Wave Modelling Project (The SWAMP 
Group, 1985) In recognition of these deficencies, the WAM (Wave Modelling) Group 
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was formed to develop a third generation ocean wave model One of the mam 
motivations of the WAM Group was to apply the mod*! latpr for FRS-i 

T he development of the WAM model involved an extensive series of investi- 
gations, including improved computations of the nonlinear transfer which governs 
the evolution of the wmdsea spectrum (Hasselmann and Hasselmann, 1981, 1985), 
the application of the improved nonlinear transfer computations to study the energy 
balance of the wave spectrum and derive improved representations of the input and 
dissipation source functions (Komen et al , 1 984), and the development of simpler 
parameterization* of the nonlinear transfer applicable m a global wave model 
(Hasselmann et al , 1985) This work was partly supported by ESA through Study 
Contract N° 6875/87/HGE-! (SC) rinally, it required the translation of these researcn 
results into a reliable numerical global wave model (Hasselmann ana Hasse'mann, 
1985 - WAM-Report). The model has been implemented at the Max-P!anck -institute 
and ECMWF, where it is used m a quasi-operational mode (Janssen et al , 1988) it has 
now been extensively tested m a senes of hmdcast cases and m quasi-operational 
forecasts. A detailed presentation of the model and its application to ten hmdcast 
stud'es is given in WAM DIG (1 988) The mode’ has meanwhile been dstributed 
together with a users* manual (Hasselmann, 1987) to more than twenty groups T 'ie 
cyc ! e 1 version of the model runs on both a CRAY XMP(ECMWF) and a CYBER 205 
(MPI) vector computer (an improved cycle 2 has recently been completed) 

Another important task which has been addressed is the interpretation and aoo" 
cation of SAR wave images The incorporation of SAR wave data in the genera 1 data 
assimilation system poses a number of problems While SEASA T provided imoress.ve 
evidence that waves can be imaged from a spaceborne SAR at 800 km height. ;he 
interpretation of SAR wave images has long been a subject of intense debate in the 
last years a reasonable consensus has emerged on the imaging mechanism as such 
(cf Section 4 1 ) The principal backscattenng mechanism in the range of incidence 
angles of the SAR on ER5- 1 is Bragg bad' scattering oy short nppies m the 5 • ’0 cm 
wave length range The SAR >s able to image ong ocean waves because me Bragg 
backscattermg by the short nppies is modulated by the onge^ waves Tne 
modulation is caused by three processes nyprodyramic interactions, c; arges n the 
local angle of incidence by the long wave tilt variations, ana sanations *n the long 
wave orbital velocity The last process gives rise to Doopier shifts n the backscatte r ed 
signal which produce azimuthal displacements of the apparent positions of the 
backscattenng elements in the SAR image plane (velocity bunching) Whereas the 
hydrodynamic and tilt modulation can be regarded as linear processes, at least to first 
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order, the velocity bunchi ng process can become strong!/ nonlinear For steep 
windsec.s tne nonlinearity can even result in a complete loss of the wave irr age 
through azimuthal smearing 

One of the mam tasks in preparing 'or the application of ERS-1 SAR wave aata 
was therefore to develop nonlinear transformation algorithms relating the SAR 
image spectrum to the surface wave spectrum Significant progress >n This question 
has been made by Bruenmg and Alpers (1985) and Bruenmg et al (19N8) The autnors 
computed the nonlinear mapping from the surface wave spectrum to the SAR image 
spectrum by direct Monte Carlo simulations Unfortunately, this method is too costly 
m computer time to be routinely applied to all ERS-1 SAR wave data collected m the 
SAR wave mode However, in the course of the present study a new closed, nonnnear 
integral expression has been derived describing the mapping of f he surface wa e 
spectrum into the SAR image spectrum. It is hoped that computations of this integral 
transform can be earned out sufficiently rapidly to be routine!) applied for the 
assimilation of F.RS-1 SAR image spectra The computation also has the advantage of 
avoiding the statistical scatter which arises in Monte Carlo simulations However, it 
s h ouid be stressed that effective methods for the inverse mapping from the SAR 
spectrum to the wave spectrum have not yet been developed Such techr,q u es w I! 
presumably involve some form of iteration, which implies that in prac r, ce SAR mage 
spectra can be usefully assimilated only if first guess wave* spectra are a /at I able from a 
wave mode! 


With the development and implementation of a reliable third genera* on g oba 1 
wave model and the clarification of the relation between the oce^n wave soectr jm 
and the SAR image spectrum, important first steps far the implementation of an 
ERS-1 wind and wave data assimilation system have now been completed However, 
further major tasks still need to be undertaken 

General overviews of the complete da'a assi mlation oroolem ha*e been gi ven n 
several position papers prepared by working groups and mdi -'idual discussants n a 
numoerof recent ESA workshops (cf Proc Conf Alpbach and Schhersee) Acomore- 
hensive summary may be found in the Global Data Assimilation Prog r amme for 
Air/Sea fluxes (an excerpt is given in the Appends) The r e ! at*on of this programme 
to the present study contract is discussed in the next section 
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1 4 The Global Data Assimilation Programme for Air/Sea Fluxes (GDAP) 


This programme has been initiated by the JSC/CCCO 1 ' Working Group on Air/Sea 
Fluxes, m response to the needs or the WOCE 21 and TOGA J) scientific community, to 
provide continuous, gndded, global fields of all physical fluxes at the air/sea interface 
(momentum, sensible and latent heat, water and solar and mfra-red radiation) The 
fluxes determine the coupling between the atmosphere and the ocean, the under- 
standing of which is an essential component of the WOCE and TOGA programmes 

The goal of the GDAP is to provide continuous 6 hourly, r x 1 global flux fields, 
beginning in the early nineties with the deployment of ERS-1, the first of the next 
generation ocean satellites In developing the programme, the JSC/CCCO Working 
Group on Air/Sea Fluxes has draw •> heavily on the contributions of other working 
groups, such as the Working Group on Numerical Experimentation (WGNE) and the 
JSC Working Group on Radiation, which are providing important contributions to 
special aspects of the overaH programme 

The programme is sub-dtvided into two phases (cf Appendix) a preparatory 
phase, from 1988- 1990, and an implementation phase, from 1991 onwards Principal 
assimilation centres are to be established in Europe, at the European Centre for 
Medium Range Weather Forecasts and in the U S , at the National Meteorological 
Center, with fu. ther assimilation systems at other national weather centres 

Within the two programme phases, fourteen major tasks are identified it should 
be noted that the present ESA study contracts address only a small sub-set of these 
tasks (projects 4 3, 4 10,4 1 2 - cf Appendix) The impact of the ESA study contracts on 
the overall programme and the inferences which may be drawn from these 
investigations for the future development of the programme are discussed briefly m 
Section 8 (Conclusions) 


n Parent bodies of the Working Group on Air-Sea e luxes are the Jo-nt Scientific Committee 
(JSC) of the World Climate Research Programme and the Committee on Climatic Changes ana 
the Ocean (CCCO) 

World Ocean Circulation Experiment 
^ Tropical Ocean/Global Atmosphere Project 
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15 Relation of present study contract to previous ESA data assimilation study contracts 

The common purposes of the previous ESA study contract, Part I (Anderson et al , 
1987) and Part II (Janssen et al , 1988), and the present contract were 

- to establish the general feasibility of using a combined wind and wave data 
assimilation system for the analysis of ERS- 1 wi nd and wave data, 

- to gather experience in the operation of such systems through experiments with 
SEASAT data, 


- to define more clearly the requirements for the implementation of a wind and 
wavedata assimilation system m time for the launch of the ERS- 1 , and 

to investigate the value of applying a wind a^d wave data assimilation system in 
support of measurement campaigns and during the ERS- 1 calibration and 
validation phase. 


Within this general framework, the individual study contracts addressed the 
following tasks: 

Pa r t I - Wind Scatterometer Data (Anderson et al., 1987) 

“ investigation o^ methods of assimilating wind scatteromete r data together wth 
the FGGE data set using the ECMWF higher resolution forecast model and 
operational assimilation system, 

validation of the SEASAT scatterometer data, 

- mtercomparison of scatterometer winds with snip wirds, 

- investigation of the impact of scatterometer winds on the analysis and forecast, 
mtercomparison of different scatterometer deahasmg algorithms 

i 
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Part M - Use of scatterometer ana altimeter data in wave modelling and assimilation 
(Janssen et a! , 1988) 

application of ECMWf- wind fields produced in Part I to a global wave hmacast, 

intercomparison of global wave hindcasts produced by wind fields *rom the 
ECMWF data assimilation system (Part I) and by w»nd fields obtamea at j PL by 
subjective analysis of SEASAT scatterometer data (Wo»ceshyn et al , 1 987), 

comparison cf global wave hmdcast with wave buoy data and SEASAT altimeter 
data, 

wave data assimilation experiments using altimeter data 

Present Contract - SAR image spectra and alt mete r wave height data assimilation 
system for ERS-1 

Most of the numerical experiments in Pa^ts I and ll of the study contract we r e 
carried out for the 1 1-day period September 6-17, 1978 Tni$ limitation was mposed 
by the expense of the atmospheric data assimilation exe'ose and the !*nrhted penod 
covered by the subjectively analyzed JPL wind data set. In contrast, most o? the 
investigations of this study contract were earned out for the full 96-day per, cd July 7 - 
October 10, 1978, for which the SEASAT scatterometer was operating Theg'obai 
wind and surface stress fielcs for this period were kindly provided by the Goddard 
Laboratory of Atmospheres. They were produced by a data assimilation exercise 
(Atlas et al., 1987) similar to that earned out by ECMWF, but using the significantly 
coarser resolution GLA modeh4° x 5°, 9 layers, as compared with 1 875 x 1 875° 19 
layers for the ECMWF model). 

The following investigations were carried out 

giobdi wave hmdeasi fur the full 96*udy 3EA3A T u^nou usmy the GlA su' face 
stress field asdnvmg field, 

mtercomppnson of assimilated wind fields and altimeter wind speeds for this 
period, 


intercomparison of the G .A surface wind and surface stress fields with the 
ECMWF and JPl fields for the common period of overlap. September 6-17, 1978, 

- intercomparison of the global wave hmdcasts for the three wind fields for the 
same common period, 

- mtercompari‘on of the 96-day global wave hindcast with SEASAT altimeter wave 
heights, wave buoys and SEASAT SAR image spectra, 

- comparison of theoretical relation between surface wave spectra and SAR wave 
image spectra with measured SEASAT wave image spectra, 

- development of a new theoretical integral transformation expression describing 
the nonlinear mapping from a surface wave spectrum into a SAR image spectrum, 

- an example of wave data assimilation for the month of August, 1978, 

- development of general methods for the assimilation of wave data in wave 
models, 

- recommendations for future actions, 

- definition of requirements for the future implementation of a general data 
assimilation system for ERS-1 wind and wave data 

The results of these investigations are presented, in the order listed, in the 

following sections. 
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Global wave hindcast for the 95-dav SEASAT period 


To test the impact of the simultaneous use of satellite wind and wave data in a wave 
model, a global wave hindcast was carried out for the full SEASAT scatterometer period 
July 7 - October 10, 1978, using as input the global surface stress field of Atlas et al (1987) 
(the wave model requires surface stresses rather than the surface winds as input) The 
stress field was constructed through the assimilation of SEASAT scatterometer wmd data 
and conventional atmospheric data using the G LA atmospheric model The global hindcast 
was carried out on the CYBER 205 at the MPI using the WAM third generation model, 
which is described m detail in WAMDlG (1988) (see also Part II of the previous ESA study 
contract, Janssen et al., 1988) ' 

In this section we compare the model hindcast with the SEASAT altimeter wave 
heights and wave buoy data A detailed comparison with SAR wave data is given later in 
Section 5, following a presentation of SAR imaging theory in Section 4 

2 1 Comparison with averaged SEASAT altimeter wave height fields 

Figures 2 1-25 show (approximately) monthly averaged fields for the GLA wind 
stress, hindcast significant wave height and direction (Custer diagram), sub-di vded 
into the contributions for the total wave field, wmdsea and swell, and the SEASAT 
altimeter wave heights (provided by JPL) The averaging periods (cf Table 1) are the 
same as the periods selected by Mognard et al (1983) (except for a data gap of 1 0 
days at the end, in October, due to a faulty tape, other gaps are due to missing 
altimeter data) 


The separation of the total wave spectrum F(f,Q', into wmdsea and swell contribu- 
tions m Figs 2 3, 2 4 was based on the criterion 


1 2- 28- u t cos i0 — 6 ) 
: -->) 



0 : windsea component 
0 : swell component 


(2 1 ) 


where c - phase velocity, 0 = wave propagation direction, 6,,, = wind direction 
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July (cf. Figs 2 1 - 2 5, Panelsa) 

11 07 1978,00:00 - 1 7 07 1 978. 16 00 
and 

24 07 1978, 10:00 - 31 07 1978, 00 00 


Auqu rt (cf Figs 2 1 - 2 5, Panels b) 

01 08 1978, 00:00 - 28 08 1978.07.00 


September (cf Figs 2. 1 - 2.5, Panels c) 

06 09 1978,00:00 - 07.09.1978, 00:00 
and 

15 09 1978, 00:00 - 30 09.1978, 00 00 


Table 1 : Averaging periods for SEA5AT altimeter wave heights 

and model hindcast wave heights shown m Fiqs 2 1-25 
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Averaged GLA surface stress fields (direction and magnitude, indicated by 
isolines and length of arrows), for August 1978 (averaging periods are given 
in Table 1) 
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Averaged GLA surface stress fields (direction and magnitude, indicated by 
isolmes and length of arrows), for September-October 1978 (averaging 
periods are given m Table 1 ) 
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direction for the hindcast swell 
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Averaged significant wave height and mean direction for the hmdcast swell 
for August 1978 (averaging periods are given in Table 1) Wave heights are 
indicated by isolmes and lengths of wave direction arrows 
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Averaged significant wave height and mean direction for the hmdcasr sw 
for September 1978 (averaging periods are given in Table 1) Wave height 
are indicated by isolmes and lengths of wave direction arrows 
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derived from the SEASAT altimeter for August 1978 
ven in Table 1 ) 








Significant wave heights derived from the SEASAT altimeter for September 
1978 (averaging periods are given m Table 1) 









The hindcast and altimeter wave heignts are seen to agree rarner wei! m the 
northern hemisphere, but the wave heights are significantly underestimated by the 
nmdcast in the high wind belts of the winter-time southern hemisphere (40'S - 60 S) 
This is seen more dearly in Fig 2 6, which shows the ratios of the hindcast to the 
altimeter wave heights. The causes of these discrepances wiii be discussed in more 
detail ;n Section 3 However, we point out here already that it appears unlike 1 '/ that 
they can be attributed to the altimeter wave height measurement o r the wave 
model, since both have been well verified under rather high wind conditions The 
problem therefore presumably lies in the surface stress field The wave heignt under- 
estimate corresponds to an underestimation of the wind stress by approximately 20 * 
30 % This example clearly demonstrates the advantage of operating a wave moae* in 
conjunction with a wind data assimilation system The wave model provides a 
valuable cross validation of two sensor systems, the scatterometer wind measurement 
and the altimeter wave height, in this case revealing an inconsistency of the 
assimilated stress field m a climatically important region of the giobe in which very 
few other independent measurements exist 

Figs 2 2-2 4 illustrate an interesting property of open ocean surface waves which 
is familiar to ocean wave researchers but is perhaps less widely known otherwse. 
most of the wave energy m the open ocean, even m the high wnd belts, is associated 
with swell, the windsea energy being concentrated locally in relatively small regions 
of high winds. This is due to the fact that swell can propagate for a few weeks over 
many thousands of km in the ocean without significant dissipation (cf Snodgrass et 
al , 1966), so that the wmdsea energy generated in rather short lived storm events s 
retained fora long time by the ocean as swe'l 

The concentration of the windsea energy m local regions is less apparent m the 
monthly averaged data »n Fig 2 3 than m the instantaneous wave he'ght distribution, 
an example of which is shown in Fig 2 7, together with the instantaneous surface 
stress field, for August 19, 1978, 00:00 G T M 

T he fact, that most of the wave energy in the open ocean is associated with sweii 
rather than wmdsea is an advantage for the assimilation of wave data m a wave 
forecast model A swell f eid, once correcteu oy a measurement update, propagates 
from then onwards m the mode! at its correc*.-d va-ue, while a wmdsea update 
r* axes back to the original incorrect value unless the wind f-eld is also properly 
corrected at the same time (see Section 5) 
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Ratio of model hmdcast significant wa^e height (cf Fig 2 2 
wave height (c? Fig 2 5) for August 1978 
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Ratio of model hmdcast significant wave height (c f Fig 2 2) to altimeter 
wave height (cf Fig 2 5) fo' September 1978 
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Fjq 2 7a Surface stress fields for August 19, 1978, 
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2.2 Statistical comparison with altimeter wave height data 


The above comparison of the averaged global fields of wave heghts was comole 
mented by a more detailed statistical intercomparison. A linear regression analysis 
was carried out between the model hindcast ana altimeter wave neignt for a set of 
8 penoas of 10 to 11 days covering the time span of the experiment (cf Table 2) 


In the standard regression analysis for a set of pairs of variables / - /. ,\ 

the regression coefficient e in the linear relation 


(2 2 ) 


is determined by minimizing the average error < v - r.r-- > w»th respect to the \ 
coordinate, where 



.V — 

i 


One obtams then the regression coefficient 


<x~> 

(We have left out the optional additive constant m (2 2), as we are interested in the 
mean factor relating the two variables ) 


One can, of course, equally well interchange x and v, minimizing the error 
< x y > with respect to the x coordinate. This yields the regression coefficient 


•> 

<v-> 
C; <xy> 


(2 4) 


The two estimates satisfy the inequality 


v 


< r 


(2 5 ) 


since 


c 

c 


< xy > ~ 

» ■) 

< X" > < V" > 


r 


1 
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wnere r represents the correlation coefficient (defined here without subtraction of 
the means, consistent with the similar definition of the regression tine slope (2 2) 
without an additive constant) 


~he regression coefficient may also be defined symmetrically by the relation 


c = u* c )- - 


: v’> 


< x* > 


(2 6 ) 


We will adopt this definition m the following The asymmetrical regression 
coefficients c c . e, may be determined from c and the (similarly symmetrica!) 
correlation coefficient r through the relations 


(2 7 ) 

c = r*r 

(2 8 ) 

r — nr 

x. 

Fig 2 8 shows the regression coefficients? and correlations r for the 8 periods 
listed m Table 2, subdivided into three regions: the northern hemisphere (2 N to 
70' N), the southern hemisphere (70° S to 2° 5) and the tropics (between 22 N and 
22° S) The abscissa axis represents the mean altimeter wave height in the relevant 
region for each ot the averaging periods 

~he statistical analysis confirms quantitative^* the general impression gamed 
already from Figs 2.2 * 2 5: the model hmdcast is reasonably consistent wth the 
altimeter wave height data m the northern hemisphere, but produces wave heights 
which are too low in the southern hemisphere by a factor of order 0 6-0 7, the bias 

. hfliAht Cm Q nn/o< an nf a twrvi/'a! craff«ar Hianram nf 

! 1 v i k. uj i > i *»>!.*> «»«>'. ■ >%. • y • "jj * • « * — ■ • r» • - — ■ »- ' T ^ — 

7/, (hmdcast) vs H< (altimeter) from which the regress.on and correlation coefficients 
were derived One data point represents an average over one day 
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CORRELATION 
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ALTIMETER - HS [M ] 
(b) 


Regression coefficients (panel (a)) and correlations (p.nel (b)) for regression 
analysis hindcast vs. altimeter wave height for the eight averaging periods 
listed in Table 2. 
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Fig, 2 9 Typical scatter diagram for //< (altimeter) v (hindcast) for one o* the eight 

averaging periods listed in Table 2 (August 1, '978, 00 00 • August 28, '978, 
00 : 00 ). 
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1) 

78071000 - 

78072000 

2) 

78072000 - 

78080100 

3) 

78080100 - 

78081 100 

4) 

78081100 - 

78082100 

5) 

78082100 - 

78090100 

6) 

78G90100 - 

78091100 

7) 

78091100 - 

78092100 

8) 

78092100 - 

78100102 


Table 2: Averaging periods used for the regression 

analysisof altimeter vs model data 
(H. and wind speed). 
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2 3 Regional comparison with wave buoys and altimeter data 


At six NO DC buoy locations in the North Pacific (4) and North Atlantic (2) 

(H Pm 1 imth rflO-M/aw .. - j — ~i . - 

' _/ “ ~ > / • • ,n,M, muuc ucuvctn uic »»ltiu flUU WdVtf 

data of the model hindcast, the altimeter and the wave buoys Tne data from three 
sources were Drought into a comparable format as follows: 


Model hindcast: The values of each model parameter (friction velocity u+, significant 
wave height mean frequency /, mean direction ?)) was determined from a 
weighted average, where the weights correspond to the inverse distance between 
the grid point and buoy location Close to land boundaries, only three or two grid 
oomts were used. 


Sdoy: From the spectral data, H s and the mean frequency f we r e recalculated 
including an H tail in the same manner as m the WAM model (cf VVAMDiG, 1988 ) 
The wind velocity was reduced to 10 m winds if necessary and then converted to u* 
values by applying the expressions given m WAMDiG ( 1988 ) in orde^ to dete^m.ne a 
suitable averaging time for the buoy data, which was available every tnree hours, the 
conation coefficients were computed for all variables using tne actual buoy values 
at model output times for noon and midnight, then the values at 3 nour$ before at, 
and 3 hours a ft er the model output times, 6 hours before, etc , up to an ave r agmg 
time of one day Plots of the correlation coefficients revealed that the opfimun 
averaging time was around 12 hours, where correlation coefficients typically peaked 
around 0 6 - 0 8 depending on the variable This result is consistent with the 
averaging incurred by the model, since the energy of a 1 0 second wave propagating 
at its group velocity would take about 12 hours to travel a distance of 3 degrees, the 
grid resolution m the model. 


Alti meter . The altimeter data set co^oca^ed wth the buoys was determined by avera- 
ging values of H s and u. over five values (- 70 km trac* 'ength) along ^ne altimeter 
track which was closest to the buoy and within a three nou r time wmaow of the buoy 
measurement A! ti meter wind soeeds were reduced to 1 0 m winds a no then conver- 
ted to tu as described above Extreme values in H < were deleted from tne data set 


Tables 3 and 4 summarize the principal statistics (the bias is the mean of the 

HlffflfOnrd tHa firrt inrl r - i - ' » • • 

- - ~ ■», *cr 9K.aK.KV> II 1UCA .->1 uit? rdlio gi me 


rms differences between model parameters and observations to the mean of the 
observations (in °/o)). In general, the agreement between the model hindcast and 
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buoy data is comparable with the agreement between the altimeter and the buoys 
The open ocean stations m the N Pacific (e g buoy 46006, cf Fig 2.12) generally 
show better agreement than the stations close' to shore (e g buoy 41001. Fig 2 ' 1), 
for which the 300 km model resolution was presumably insufficient A comparison of 

. i . . - _ _ - T. 1 11 1 11 c. rtnoitc that 

the triction velocity ana wave iieiyiiuimocnei ■ 1 «• «- - — 

much of the discrepancy between the different instruments can oe attributed to 
local, relatively short-lived wind events whirh. nr rhe one hand, are not adequately 
resolved by the altimeter and model and which, on the other hand, are not 
representatively sampled by a point measurement Although the general consistency 
of the intercomparison gives support to all three measurement systems, the scatter is 
nevertheless greater than is typically found in model-measurement mtercomparisons 
involving either larger averaging areas (cf Section 3, for example) or point 
measurements accompanied by a higher resolution determination of the 
surrounding wind field (cf. hurricane and North Sea hindcast studies discussed m 
WAMD'G, 1988) 

From the viewpoint of an optimal measurement strategy for the validation and 
calibration of ERS-1 during the initial commissionary phase, as well as the ong term 
monitoring of the satellite sensor system, we therefore recommend giving higher 
pnonty to (i) global verification methods, based on combined scatterometer, 
altimeter, and SAR wave mode data, as outlined in Sections 2 1,22,3 anq 5, and 
(ii) measurement campaigns, in which in situ measurement s’abons are deployed m 
arrays, rather than as single stations. However, where available, long term Single 
station buoy measurements can nevertheless usefully augment these othe r 


approaches 
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TiWH Statistics ot model vs buoy intercomparisons for individual buoys (to be continued next page) 





BUOY 46001 BUOY 46005 BUOY 46006 

56 N 212' E 46“N 229"E 41°N222 E 

Variable (56 observations) (145 observations) (145 observations) 



**3 


Tab[e_4 Statistics of model vs buoy intercomparisons for individual buoys (continued from last page) 
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Positions of the six NOOC buoys used for model hmdcasl, altimeter and buoy 
data m ercompar son 



288 00 



JULY 1978 


AUGUST 1978 


SEPTEMBER 1978 


Tim* series wind and wave data from wave hmdcast altimeter and wave buoy 
at location of buoy 41001 in N Atlantic (cf Fig 2 10) 
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BUOY 46006 (41 00 N, 222 00 E 


ALTIMETER o 
MODEL x 
BUOY 


FRICTION! velocity direction 




JULY 1978 


AUGUST 1978 


SEPTEMBER 1978 


■Tme ser.e* wind and wav* data from wave hmdcait. altimeter and wave 
buoy at location of buoy 46006 in N. Pacific (cf. Fig 2 10). 
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3 mtercompanson of wind he'ds and wave hindcaV 's 


3 1 intercom panton of Gl_A surface wmd and stress fields with ECMWF and JPLA^ /jare 
wmd fields 

To determine the possible origins of the discrepancy between 'he model hmdcast 
wave heights and SEASAT a*timeter wave heights m the sourthern nemisohere. an 
mtercompanson was made between the GLA surface wind and stress fields and two 
other wind fields which were availaole for sub-periods of the 96-day hmdcast period ■ 
the assimilated ECMWF 1 000 mb wind fields for the period September 6*17 (^nderson 
et al . T 987), and the subjectively analyzed JPL 19 5 m height wind fields for tne penod 
September 7-20 (Wo»ceshyn et al , 1987) Details of the wmd and stress data sets o *• 
given m Table 5 

in the rol lowing two sub-sections we consider first the mtercompanson tor the 
three wmd field data sets and rhen the GLA stress wmd f -e'd mterconoar-son it was 
found that the mam problem appears to lie m the GLA southern hem.sone^e stress 
held analysis, while the three wmd field products are reasonably conyste v fo' ai- 
regmns of the globe This finding is borne out by the mtercompar son o' tne "hree 
wave hmdcasts for the period of common overlap discussed m Section 3 3 

All following mtercompansons are expressed m terms of the symmetrised 
regression coefficient c relating pairs of fields U. f v » in accordance with equ (2 6 ) 

3 l 1 Wmd field mtercompansons 


Table 6 summarizes the results for the three wind field regressions for different 
pairsof the GLA. ECMWF and JPi wmd fie l ds m adamon to me gioba 1 stansncs data 
are also given for the northern and soutnern Hemispheres ana me ‘ r oon.$ ^ gs 3 1 - 3 3 
show the dependence of ? and r on wma soeed in each of the three ia f tuc r at be its 

* .9*- m **«**<’» . n * 

; or eacn of the three mtercom^ * w acr 'mr m me se gures vor'esoo* i us .0 an 

average formed for *he ccrresponc • - a oart cuia f a n a ys s me 

Trie principal concision to be drawn from these r esuts s that tne three wmd 
fields show a surpnsmg level of agreement the regression coefficient lies dose to one 
for all three mtercompansons m al! three regions, and *he correlation coefficient is of 
the order of 90 - 95 % in all cases it appears that the discrepances between the 
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aitimeter wave heights and the moael wave hmdcast m the southern hemisphere 
shown in Section 2 1 cannot be attributed to the wind field analysis The problem lies, 
in fa^t, in the southern hemisphere GLA stress fields, as shown ,n the following sub- 
section 

a similar level of agreement is also found m the wind directions As example, 

Fig 3 4a shows the global histogram for the difference m wind direction AO between 
tne GLA and ECMWF wmd vectors Most of the errors are less than 20 : Fig 3 4b shows 
the distribution of the mean wind speed and the standard deviation of the wind speed 
as runction of the directional deviation As to be expected, ! arger directional 
deviations are gener - !ly associated with weaker wind speeds Similar results are found 
for the other two intercomparisons and for the breakdown into the three latitude 
zones 

3 1.2 Intercomparison of GLA surface wind and surface stress fields 

To test the consistency of the GLA ( 1 000 mb) Surface wmd speeas . With the GLA 
surface stress values i, a regression analyse was earned out be'ween and dp lt whe r e 
o , denotes the dens.ty of air Thp rpgrpssmn rnpffioenr ” then represents an estimate 
of the average drag coefficient cp Since the local (temperature dependent) p . values 
were not available on the GLA data tape, p a was taken as constant^, - 1 23 kg m3) 

A problem arose in the assignment of time levels. The wind field analysis at time i. 
represents an instantaneous field obtained by combining the first guess wmd field 
provided by the model forecast, derived from the previous analysis time /V/, 6 hou r s 
earlier, with conventional atmospheric observations and scatterometer winds at (or 
close to) the analysis time t„. The surface stress fields, on the other hand, are computed 
as the average vaiues of the model over the period r.../ to As perhaps to be 
expected, the best correlations were found between u t.j p, : and o 1 ' t.. t . This is 
illustrated in Fig 3 5, in which the giobal regression coefficient rand correlation r is 
shown as a function of the lag Ar = 1 . 1 r ,p„ ) - r > -) The correlations discusseo in the 
Allowing were made for At = 6 hours 
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(M/S) 


■ 60 . 

THETA(DEG) 

(a) 


- 180 . - 120 . 


■ 60 . 0 . 60 . 120 . 180 . 

THETA(DEG) 

(b) 


Probability distribution of devation in direction AO between wind vectors 
(pane' (a)) and mean wind speed and rms wind spread (panel (b)) as function 
of AO for GlA-ECMWF mtercompanson (70 : S - 70^ N) 




The cause of the underestimate of the southern hemisphere wave hmdcast is 
revealed m Figs 3 6a, b These show the dependence of the regression coefficient c 
(j e. the drag coefficient) and the correlation coefficient on the wind speed for the 
northern, southern and tropical regions. We expect an increase of r with wind soeea .n 
accordance with Charnock (1955) or other similar drag law relations (e g Wu, 1 982) 
instead, "decreases with increasing wind speed The northern hemisohere values of c 
of the order of 1.2- 1 7*10 ^ are reasonably consistent with emomcai drag laws for 
the moderate wind speeds characteristic of the northern hemisphere summer T n e 
values for comparable wind speeds in the tropics are generally too low However, the 
biggest bias is found in the southern hemisphere For mean southern hemisohere w nd 
speeds of the order of 1 1 - 1 2 m/s, Wu's formula gives co = ( 1 52 - 1 58) • 1 0 3 which is 
20 % to 60 % larger than the "values in the range 1.0-1 2- IQ' 3 For a fully develooed 
spectrum, the significant wave height in the WAM model is proportional to the surtace 
stress. Thus, it appears that the bias found m the model wave height hmdcast relative 
to the altimeter wave heights in the southern hemisphere, which was of about the 
same magnitude, can be explained by the underestimate of the GlA surface stress 
f elds m this region 

We can offer no simple explanation here for the cause of this underestimation 
However, it appears probable, in accordance with a suggestion of the GLA group, that 
»t is related f o the 6 hourly averaging involved in computing the stress field Through 
the nonlinearity in the definition of c, equ. (2 6), this leads to an underestimation of 
the drag coefficient, and a similar nonlinearity of the wave model response leads to a 
comparable underestimation of the wave height The mam conclusion we wish to 
draw from this analysis is that the operation of a wave model in conjunction with 
general data validation cross checking procedures, applied to both model products 
and observed data, can rapidly identify and locate problems m different data sets The 
operational application of such procedures would clearly be ve r v valuable for 
continually monitoring the performance of the ERS-! wind and wave sensor system 
and algorithms 


Other problems related to the original scatterometer wmds (which we^e nor 
considered in the present project) have been identified m the orevious ESA stjdv 
contract, Part I (Anderson etai , 1987) and Part il (Janssen et al , *988) AH three 
investigations demonstrate that the simultaneous application of sopntst.C3 f ed mode s 
a general data assimilation system and a wide spectrum of data validation techniques 
on an operational basis is an essential pre-requisite for the implementation of ?■ 
reliable end-to-end system for ERS- 1 wind and wave data processing 
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3 2 intercomparison of GLA surface wind field with SEASAT altimeter wind speeds 


For completeness we show in this section an intercomparison of the GLA surface 
winds and the SEASAT altimeter wind speeds, which were obtained simultaneously 
with the altimeter wave height data shown in Section 2 1 The format of the analysis is 
identical to the mtercomparison of the altimeter and wave model hmdcast wave 
heights presented in Sections 2 1 and 2 2 


None of the analysed wind fields considered m this study actually made use of the 
altimeter wind speeds It is generally accepted that the al nmptpr w/mrt which 

are inferred from the altimeter measurement of the rms slope of the sea surface are 
ess reliable than the scatterometer wind speeds derived from the short 8ragg 
scattering surface wave ripples This is confirmed m our analysis 

Fig 3 7 shows the global distributions of time averaged altimeter wind speeds for 
the three approximately monthly averaging periods listed in Table 1 This may be 
compared with the corresponding F.g 3 8 for the GLA wind fields The altimeter winds 
are clearly significantly lower than the GLA winds, particularly m the high wind regions 
in the southern hemisphere Since it was shown that the GLA wind fields, in contrast to 
the GLA stress fields, are consistent with other data m both hemispheres, the 
discrepancies must be attributed to the altimeter winds rather than the GLA wind 
field 

An underestimation of the w>nd speed m the SEASAT altimeter algorithm at high 
wind speeds is found also m the statistical mtercomparison of the GLA and altimeter 
wind speeds shown in Fig 3 9 The high wind speed bias is of the order of 25 % A 
recalibration of the altimeter wind speed algorithm (which should be directly 
applicable also to the ERS-1 altimeter) appears called for Smce the altimeter measures 
therms sea surface slope, which is known to be a function not only of the local wind 
speed, but also of the sea state, an improved algorithm should preferably combine 
wind and wave information 


58 



averaging 
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Of August 1978 
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Averaged global field* fv. SEASA7 altimeter wind speeds for September 1978 
(exact averaging period* are given m Table 1) 
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■SL. Mi Averaged global fields for GIA 1000 mb wind speeds for Sep 









3 3 intercompanson of wave node! hindcasts for GLA, ECMW C and jPl surface 


f orcing 

it is of interest to investigate how tne comparison between the surface wind and 
surface stress fields for the three wind field analyses discussed above is translated into tre 
corresoondmg wave model hindcasts As has already been pointed out, wave model 
simulations gene r al I y provide a sensitive test of the input w md fields 

Rather than reoeat the statistical analysis presented <n Section 3 1 for the wind soeecJ 
and directions now in terms of the corresponding wave variables, the significant wave 
heigpt and mean propagation direction, we show in Figs. 3. 10 * 3 12 a comparison of the 
three global hmdcast fields for the significant wave height, averaged over tne common 
period September 6-17, 1978, of the three hindcasts. 

^he ECMWF and JPL wmd fields were transformed .nto stress fields requ»rea as input 
f or the wave model using Wu's formula (which was also used m the cal'bration of the wave 
model, cf Komen, 1985. WAM-QiG, 1988). 


'• = c nP., ( -' 


wre r e 


r, 


1 2*75 10-!, 


‘0 *+ 0 065 U> l0-i- 


L‘ < 7 5 m s 
i 5 : 7 5m s 


No distinction was made between the 1000 mb ECMWF winds and the 19 5 m height 
JPL winds. 


p gs 3 10-3.12 support the general conc!us<ons of Section 3 1 tnere $ r easonaoiy 
good agreement between tne three forcing fields m the nortnem nemiso nere, while tne 
GlA forcing is clearly significantly weaver than the other two forcing fields m the mgh 
wind region of the southern nemusohe r e A comoanson with tne altimeter wave heights 
shown m Section 2 i for the entire 96-aay pe' od suggests that the ECMWF and jPl forcing 
£ ieids are probably of about the right magnitude m the m.d- ana hign-latitude soutnem 
hemisphere, while the GLA stress field is too weak in this region 


Apart from these general relations, the figures reveal a number of regional 
deviations, which will not be pursued here, but would clearly be of interest m applying a 
wave model for wind field validation in an operational setting 
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Significant wave height for -wave hmdcast for JPL surface wind field for the 
period September 7-17, 1978 







Bass et a\ (1968) and has been extensively tested and verified m field and wave 'ark 
experiments On the basis of this model, a rather complete tneory of SAR imaging of a 
random ocean wave field can be developed (cf Aipers et ai . 1 9 79, Kasselmanr et a! , 
1985) 

An important feature of this theory is that SAR imaging can oecome strongly 
nonlinear Although the hydrodynamic and tih modulation can normally be regarded 
as linear processes, the "'nr, on effects become strongly nonlinear *n typical w ndsea 
situations The doppler shift arising c rom the radial component of the orbital velocity 
of a b<Kkscattering facet produces an azimuthal displacement of *he apparent position 
of the backscattenng facet m the image plane. For displacements small compared w»th 
-he characteristic wavelength of the long waves, this ‘velocity bunching’ mechanism 
can be treated as a linear process It can be characterized by a ve'octty bunching 
modulation transfer function wh<ch, just as the hydrodynamic and hit mtfs, 
contributes to the wave image contrast. However, if the displacements ~o:ome 
comparable with or larger than a quarter wavelength, the ve'octy bunch*rg 
mechanism becomes nonlinear, and for still larger displacements ultimately resets n a 
complete smearing of the image 

'he transformation or a surface wave spectrum into a SAR ,mjge soect'um ras 
oeen computed numerically in the weakly and strongly nonlinear regmes by Monte 
Carlo techniques (Aipers and Brunmg, 1986, Brunmg, 1987, Brunmg etai , 1988) From 
an extensive compilation of such simulations for afferent wave conditions and SAR 
parameters. Brunmg et ai (1983) denved general relations for the distortion of the 
surface wave spectrum in the SAR image spectrum To first order, the principal effects 
are that short waves travelling m the azimuthal direction are filtered out, the soectrum 
sskewed towards the range direction, and the mean wavelength is increased 

~he Monte Cano computations are based on individual reaiizahons of a r^naom 
surface wave f e } d c o r eac H realization, the .ndiviauai movmg oac- scatte r rg surface 
elements m each scene are maooed mto their corresoonamq positions m me SAR 
■mage plane, ana the resuming image s Founer transformed "he SAR .mage variance 
spectrum is then computed by ave r agmg the spectra for an ense^ote of suo 
realizations This orute : orce metnoa ■$ rather expens, ve m computer tim* 
cannot be applied operationally *n a wa*e model to rebate the model output soeevum 
to the large number of SAR image spectra which will be provided continuously by 
future oceanographic satellites such as ERS* 1 or RAOARSAT 



In this section an alternative computation is proposed it is shown that the wave 
spectrum and SAR image spectrum can be related directly m both the linear 0 nd 
nonlinear mapping regimes, through a general nonlinear integral equation The 
transformation can be readily evaluated numerically and can therefore be used in 
routine operational applications 

At first sight it appears rather surprising - in view ot the notorious closure p r obiems 
of strongiy nonlinear systems * that a closed expression can be derived relating the 
surface wave spectrum and the SAR image spectrum in the strongly nonlinear reg me 
However, closure is possible in the present case because the input and output fields do 
not interact dynamically, and to a good first approximation the mout surface wave 
field itself can be regarded as linear and Gaussian 

Similar techniques to those developed in Sections 4 3. 4 4 have been applied 
previously in the computation of the dissipation of finite-depth surface waves by 
bottom friction (Hasseimann and Collins, 1 968) and m the determination of the 
structure of the dissipation source function due to white capping (Hasseimann, 1974) 
The basic concept is to consider the isolated impact on the process of the infinitesima 1 
contribution from a single Fourier component of the input field This can be 
determined rigorously, because for a Gaussian wave field the infinitesimal 
contribution to the continuous spectrum are statistically independent The resulting 
(nonlinear) mtegro-differential expression relating the infinitesimal change ,n the 
surface wave spectrum to the associated infinitesimal change m the mage spectrum 
can then be integrated analytically to ye'd the general nonlinear integral 
transformation from the ocean wave spectrum to the SAR image spectrum 

In the following Section 4 2 we summarize first the relations for RAR (real aperture 
radar) and SAR imaging, as derived from the standard two-scale scattering theory and 
introduce notations The general nonlinear transformation relating the SAR image 
spectrum to the ocean wave spectrum is then derived m Section 4 3 A mo'e detailed 
discussion of the various terms occurring in the transformation is t, ven m Section 4 4 

4 2 Mapping of the ocean surface into the SAR image plane 

The SAR surface wave imaging process may be divided into two processes the RAR 
(real aperture radar) imaging processes, consisting of hydrodynamic and tilt 
modulations, and the motion effects, which are spec fic to the SAR and have no 
influence on RAR imaging 
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vVe consider first the RAR .magmg process 


Decor" cos' ''g the surface «!*vaton i- x, t). ore the loco car: i,:a*fe..!'.y ny» 
section oi x, t sensed by a RAR .nto 0 scre*e Four-er senes, 
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(4 2) 


vs here o os the soatiaMy averaged mean (specific) cross section and ^ the cross-sect. on 
modulation factor, the standard two-scaie scattering rneorv for hyorodynamic and rot 
modulat'On yields a 'inear re'afion 


'o = r k\ ,J » 

between the Fourier components of and with a net RAR comoie* modulat on 
transfer function (MTF) T k - T k h - J - Tj which consists of me sum ot the 
hydrodynam*c and tilt modulation transfer functions For our purpose* we rega r d l 
as given (It is notationally more convenient to work with discrete Fourier sums than 
with continuous integrals as we snail be considering later differentials with resoer to 
maividuai Fourier components he discrete representation avoids the rather 
cumbersome functional derivative notation ) 

in equ (4 1) the frequency isgiven by the free gravity wave dispersion relation 
u/ r ■*, where j is the acceleration of gravity it should be no r ed that equations 
(4 1), (4 2) represent three-dimensional w. /enumoer-freauency soectrj wnicn have 
been reduced o tw o-d-mensional soectra because me rr equenoes are co n st r aned to 
i'e on thea'soe's.o'* 1 solace m contras: to the fwo-c^e^.o^ai wavenumber 
spectrum of a 'rozen surface cons.dered be'ow (equ (4 4)) =our<er comoonent* 0 r 
Cf.posite sign ip k represent waves *rave!l»nq m oooosite directions ana are not 
re'ated 


Radars, nn the omer hanc. produce quasi-mstantaneo^s images or me surrace. so 
that the relevant Four-er decomposition of the image i r at a fixed time, t = if. say 
takes the form 


(4 4) 


wnere 


/' ! ) - I l - ^ / >>xt) i fk 
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r 1 1 


and the image intensity is normalized with respec' to the mean mag- intensity (No'e 
that the explicit complex conjugate term of eqns (4 1), (4 2) is m vmg here, as it s 
already included in the sum over positive and negati >e k ) 


p ° f a RAR, the 'mage intensity »s simply prooortional to 'he c"Oss section, so that 


/ i r 


1 r. i/'.j 


(4 6 ) 


ana equa^ons(4 1) - (4 5) yield 


F'om equ (4 7) we find then that the relation between ‘he ocean wave and mage 
variance spectra F iy P k , respectively, defined by 


V>.V fk . s N.^ • 
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where *ne cor^e r ec oa^enmeses < 


> delete e^semote a . "ages, i$ g ,en ov 


4 k 



(4 


’ h * RAR image / ' r corresooncs to 'he .mage wn.cn , $ar would oroduce . ' «ne 
tea surface were frozen Consider now tne distortion •„>* this image induced by tne 
motion effects T hese consist of two ‘e'.-is: 30 armuthal d.spiacemen' 5 , of tne 



position or a backscattenng element in the image p‘ane, and an azimuthal smear, ng or 
broadening h.i of tne (theoretically infinitesimal) backscattenng element m the image 

“he azimuthal displacement is proportional to the range component r of the 
ave r age au ^action (i.e orbital) ve'octy of the back scatter, ng e ement (cf Aipe^s ana 

R:j tPH^rh^ 1Q7Q) 


£ = lit 


KP) 


where 


and p '$ the slant range, i the velocity of the SAR The orbttai v/eioc*tv »_• •$ oe fir 'eo 
as the time average velocity over tne period during whicn tne scattering e'ement is 
v-ewed by the SAR Normally, this is small compared with the wave oer.oa, so that. f o 
first order, we may set c equal to the instantaneous orbital ve'ocfy component m the 
range direction at the median time assigned to the SAR ,mage 


According toe andard surfar va e theory (cf Philips, 1977), 


r = ^ T^ »xp'.i\kr- (o/|) cr 
k 

wnere tne range velocity transfer function is giver by 


(4 13) 


T = < — trn.sf) — d k PkO. 


(4 U) 


and i), h v aenote the radar incidence angle and r he wavenumoer comoonent m tne 
horizonta 1 racar too 1 : direct'd", respectively 


“he smearing term <S.v ^ is normally represented as the sum of an acce ; e r ation :e r m 
'•x " and a ve'ooty soreac + e rnr '\v < c t Hasse^ann et ai . 1 985 ) 

^he acc9'e r at : on :e r m . is a second order term ar s*rg from rhe variation of tne 
instantaneous oroital velocity component ;• during the SAR viewing interval This yteids 
slightly different effective displacements £ for the beginning ana the end of the 
viewing period The tern is generally smaller" than the velocity spread term (cf. Aipers 


and Rufenach, 1979, Alpers et al , 1981, Hasseimann et al , 1 985) and, to simplify the 
analysis, will not oe considered m tne following 

^he velocity spread term *'..r ^ is icentirai m physical origin to the azimuthal 
aisDlacement term it 'S distinguished in SAR theory from tne disoiacement term on ! y 
'or formal reasons. In a comomed microwave backscattermg and SAR imaging theory, 
it is convenient to distinguish between c/ o separate two-scale models, the 
hydrodynamic/backscattenng and the SAR two-scale model (cf Hassse'mann e f ji 
1 986) In the hydrcdyn lamto'backscdttef ing mudei, which we have been consider! ng so 
far, the separation scale L^. {4 corresponds to the dimensions of the backscattermg 
facets This must be large compared with the wavelength of the Bragg scattering 
waves but small compared wth the wavelength of the modulating long wave fieid 
Typtcaiiy, v is of order 1 m. in the SAR two-scale model, on the other hand the 
separation scale L ^ 0 i s aefi nea as the SAR resolution scale, wmch is tyoicaily of oraer 
29 m The SAR ts unaoie to distinguish between individual backscartermg f acets wtnm 
a SAR resolution cell, and therefore maps the entire enserroie of back scattering facets 
w'thm a resolution ceil into a smgie image pixel This is azimuthaUy displaced by an 
amount At = [k determined by the mean orbital velocity of the facet ensemble r he 
devationsox = L - £of the individual facet displacements from the mean value "for 
the resolution cell then results in a smearing of the image of the resoiuhon re 1 1 Tne 
rms value of this smearing is given by 



~he term is essentially determined oy the contribution to the rms orbital velocty 
from the region of the wave soectrum lying between the two separation scales /. 

and V? 

in the oresc ontext it is moortant to note that me introduction of a velocity 
spread term is a *ormai.tv. wmcn m some cases nas conreotuai advantages, out is not a 
necessity. In our case it is more convenient to work entirely n the framework of the 
basic nyd rodynamic back sea tte r: ng two-tea) e mode) , w« thou t i n vo* i ng t. e vei oc tv 
oreaa concept. m this cc^jre .\e rh^ > mace as composed of a 

superoosi tion of the independent images or individual oas*c backscattermg facets, 
rather than of individual resolution ceils. 

Each backscattermg facet experiences an azimuthal displacement m the image m 
accordance with (411) ar.d some defocussmg through the face' acceie r ahon (which 


we ignore) in addition, the individual facet images will then be smeared by the finite 
resolution of the SAR. Thus we consider a pure velocity bunching' theory with explicit 
velocity spreading' To simplify the analysis, however, we ignore also the smearing oy 
the finite SAR resolution and assume the SAR has infinite resolution The effect of the 
finite SAR resolution can be easily taken into consideration, if necessary, as an 
additional filter at the end of the analysis Qualitatively, the acceleration smearing 
term can be similarly handled, although its rigorous inclusion m the theory is 
somewhat more complicated The effects we have retained do, m fact, represent the 
critical processes limiting the ocean wave imaging performance of a SAR m the 
nonlinear imaging regime 

Applying equ. (4.1 1) to a continuum of facets, we obtain then as the relation 
between the SAR image and the RAR image in the present velocity bunching model 

C ir) = | / R < r'lbir - r' - f.'r’n dr' (4 l6) 

where £ = af. and a denotes the unit vector m the azimuthal direction 
Integrating over the 5- function, equ (4 16) yields 


o or 

Mf) dr K = r -W> 


where the Jakobian 


■ dr' i c^lr’l (4 18) 

! — ; = i + — - 
i dr i or 

The 'velocity bunching factor’ ! dr * dr 1 represents the va nation m the elective 
density of backscattering elements in the image plane resulting from the compression 
or dilatation of the originally homogeneous distribution of facets through the spatial 
variations of the facet azimuthal displacement. This enables the SAR to image ocean 
waves even in the hypothetical si tuation in which the RAR transfer function vanishes 
i.e. ( H <n = const 

For the case 
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4/r 

the velocity bunching factor can be expanded ,n a geometrical senes and truncated 
after the imear term if it is furthermore assumed that the modulation of o»s small, as 
required for a linear PAP imaging theory, the dom.nant term m {4 t 7) «s tne nr.ean zeo 
velocity bunching f actor, anc one ootams tor tne SAP image amplitude spectrum, m 
accordance with (4 4), (4 13) 
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where * n e velocity ounchmg modulation transfer function 
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hus m the linear aporoximation 
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(4 22) 


an c me mage variance spectrum is given by 
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k 



where the SAP imaging transfer function 


(4 23) 


-tO r>t? -r> ') ( ** ^ 

k ~ k ^ k 

'he condition {A 19) is generally satisfied for swell spectra However, in many 
situations, for example for short windseas, the inequality does not hold or is even 
reversed (cf Hasseimann er al , 1985, Bnjnmgeta! 1988) intmscase, equations 
(4 1 7), (4 13) represent a strongly non.: near transtO" mahon and we must seek 
alternative method? for deriving the relation between tne surface wave and SAP 
image spectrum 


4 3 The nonlinear mapping from the surface wave to the S/-* R image spectrum 

he relation describing the nonlinear transformation from tne surface wave 
spectra to the SAR image spectra is derived m two steps. First a differential relation 
describing the effect of an infinitesimal perturbation in tne surface wave soectrum on 
the SAR image spectrum is deveiooed. This expression is then integrated to yield the 
final nonlinear transformation relation 


he Fourier component (^> of the SAR image :s oo tamed by tak mg the Fourier 
transform of (4 17) 


/?. 1 | : - 

k A ' .lr 


•'xp < — :k-ri= — | i/r'/irVtpl - ;k(r’ - £,',r'/)| 


A 25) 


Here A denotes the finite area of tne sea surface corresponding to our discrete founer 
reresentation (in the final result, we shall, of course, let A — *). 


Substituting the Fourier representation (4 4), (4.7) for //? in (4 25), this becomes 
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(4 26) 


We obtain then for the variation o/^mducea oy perturbing tne set of Fourier 
comoonenrs of themput wave heid by the .nfmitesimal increments .it • 
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where 
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The summation term in the second parentheses arises from the differentiation of tne 
exponential factor *>rn( -ikfy r r > I in (4 26) 

Consider now the perturbation of the variance spectrum F^ arising from the 
perturbation (4.27). At this point we make the essential assumption, which will enable 
us to close the moment equations, that the perturbations cic*, for different wave- 
numbers are statistically ^dependent of each other and are aiso statistically indepen- 
dent of the existing wave field. The assumption is justified by the Gaussian property of 
a surface wave field, which may be regarded as a superposition of an infinite number 
of statistically independent infinitesimal wave components. We shall follow up this 
concept explicitly later by regarding the surface wave field as being graduall> built up 
from a zero sea state through the continual addition of statistically independent 
spectral increments. 

Under the stated statistical assumptions, we obtain immediately the differential 
relation: 



k* 


3 e*ore turning to the integration of equ (4 29), we must first evaluate < 

The essential nonimear terms in this computation arise from exponential expressions 
of tne form o xp ■ $ir w >) > Since £ is a linear functional of c. its probability 

distribution is Gaussian, and the expectation value of the exponentials can therefore 
oe determined analytically. The computation of > given in (4 5) 

(Appendix). 


Let us now construct the prescribed surface wave field by building uo the spectrum 
F(k) from zero to tne final state through a linearly increasing continuum of wave f.elas 

A A 

F k \, = \F k wnere 0 - \ - l,andoF k -F k o\. Denoting the associated SAR image 

A ^ A 

spectrum ana matrix .W kk , in the intermediate state \ s.rnttan y as F k '(\), V/ kk » i \), 
respectively, equ. (4 29) then becomes 
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(4 30 ) 


This can immediately be integrated with aspect to V to yield 


(4.31) 


F s = V 

k — 
k' 




where the transformation matrix 


T = 
kk* 


< 
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The computation of k - is also given in the Appendix. 


(4 32) 


The final mapping relation (4.31) between the RAR and SAR image spectra is seen 
to be a quasi-linear transformation. The transformation matrix T kk . is in general a 
nonlinear function of the wave spectrum. In the limit of a very small wave spectrum, 
the relation becomes linear: T kk . becomes independent of the wave spectrum and 
diagonal. As the wave spectrum is gradually increased, T kk . develops off-diagonal 
terms and the matrix elements with high azimuthal wavenumber components are 
attenuated (azimuthal high wavenumber cut-off) This is discussed in more detail in 
the following section. 


4.4 Structure of the transformation matrix r hlt . 


In Section 4.5 ( Appendix) it is shown that the transformation matrix consists of the 
sum of two terms, 


T ~T' + T n 

kk kk' kk’ 


namely the 'filter transformation matrix' 


(4 33) 
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(4 34) 


and the background transformation matnx' 


k,k' 


= 8 k.K'-k ,P,r,+ 7 ' l <r,)+C k.k-. 


•r.,(r) 


(4 35) 


where the coefficient matrices .4 k k . v 8^ k . k and C’ k k< k depend on integral 
properties of the wave field (cf eqs (A 20) - (A 22) 
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The structure of the transformation matrices becomes clearer when it is recognized 
that the functions ^(y), ^(yjand <My) eqs (A. 15) -(A. 17)) act as high azimuthal 
wavenumber cut-off filler functions which control the form of the coefficient matrices 

■Wk and fi k.k-.k- 

consider first the h'lter transformation matrix r f kk , Fory-0, / and equ 

(A. 11) yields A k k , k — o k k „ Thus in the limit of a very smail azimuthal wavenumoer h x 
or very smail orbital velocities, this term reduces to the linear SAR imaging expression 
(4.23) As y is increased, the function .-\ k k * k is broadened with respect to the second 
variable k* - k. The wa^e component k is thereby spread in the image spectrum into a 
band around k, and the energy of the baoj is attenuated. The spreading and 
attenuation increases with increasing azimuthal wavenumber k x and mc r easmg 
orbital velocity. 

The second term, the background transformation matrix r^., vanishes entirely as 
the orbital velocity approaches zero T 'usm this limit equ (4 33) reduces to the RAR 
imaging relation, as, of course, it rv :t. For finite orbital velocities, the term T° k k> 
yields a background continuum For \ — 0, «J>Jy) - 1 and B. ... ^ , <\ - 

so that the continuum degenerates m this limit to tne quadratic sum and 
difference wavenumoers of the wave field, as expected for the lowest order auad'atic 
interaction of a nonlinear process For finite the function B k k . R is broadened with 
respect to the second index variable, so that all wave components contribute to r ne 
continuum. 

4 5 (Appendix) Computation of <UVf. . ,j*> and 7\ . 

— ■ i — -i ■■■■» — ■ . nh i i ■ JuL 

From equ. (4 28) we obtain: 


‘^kk’ • ” " drMr’Vrph ik*- k< - /k rT — . 


i r k 1 r k . k-r*> 

k” ‘ 


(A 1) 


-tVj r“.j ^ /*• apuikY-kVn 


kT k” 


where P ? k as m ecu -7) 

To determine the exoectation value <|Af kk .|^> we expand the exponential factor 
t’xp i -tk rfc r f ) - $( r w> h with respect to a given pair of infinitesimal Fourier amplitudes 
^ k m ana their complex conjugate amplitudes c*_ k «, c* k «*. vVe nave nrst 

Or') — ^vr M "i = A£ = uk”r') - apukVTlf >•> ex pi uk m rh - e.vp wk"V*>* 

(A.2) 

■+■ &T U k „C expt/kVi - expwk’V'h) +■ c _ kt Je.Tpwk , V> — *-rp<ik"V)i-’- Ay ^ 


= Ac> Ac 


'•t'.sr 


where Ac contains the four infinitesimal Fourier components of interest and Acres* * 
the remaining field which contains no components at these wavenumbers 

Thus 


to i -ik AO = expi — rAr AO > -<1 - A£ - tf“< 2 iA£ )' - . 

' r*st r r 

The cor-elation of this expression with the amplitudes *, k «, c k - w\ k ». c\ k ~ occurring in 
(A t) then yields 


.Vf kk . j > = A -1 j r/r raplnk’ - kir! 


(A 4) 
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where 


r iri = rf * P iri 
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1 s inaeoendent of t ne exoansion of the exoonent 
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( r-r* - r". r" = r () ) 
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is derived from the linear term in the expansion of the exponential function, 


r <r> — — T' -PiD-Pl-ri 
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is derived from theauadratic terms in the expansion of the exponential function, 
and 
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(P^ k as in equ (4 10 ) 
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in deriving (A. 4) * (A. 8) we have set r = r* r" in (A. 1 ) and carried out one integration 
over r‘ for fixed r, noting that for a statistically homogeneous wave field the integrand 
depends only on r. We have also replaced Oxp * by <*xp - <kA£ >, since 

A (>^*1 ana At deviate by an infinitesimal quantity 

\f w k . $ obta-ned from Af kk .by replacing £ k . c k by \ T - c k and \ y * £ k , respectively 


>Vf \ > ! > = A~ l dr *xp\i\W - k srh " — ;k\ A£<r? > 


i ! r; S 1 - \-<P(r) + T »r»<K\-T <n 


(A 9) 


We need to determine still the expectation value of the exponential function, 


r, m = < rxp — r k.\ r ‘ A£j > = V.vp - ( k \-q > 


(A 10) 


wnere 


q ~ i.Ut-u-’) - t’»r M )| 

The variable q is Gaussian with zero mean and known variance 
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: n’>=r]Li T u lF k~ ' ’ r Lk! F _»<f l ~ roskr 


so that the function <f> can be readily evaluated : 

f > (A 1 

<f)ik. r, \) = *\r p i — «’\ < rp >2 ) 

All terms occurring in equ (A. 1 ) have now been expressed directly i r, terms of the 
surface wave spectrum 

To determine the matrix f kR . = / <!.Vf kk . I- > d.\ occurring in the final expression 
(4 31) we need now only to integrate equ. (A. 9) witr respect to V This yields 


= A “ ) dr ^xp l (k‘ — k) ■ r> • 4> p k, r J T k 


<k.rh (P< r) ■*- r^irM <{> ^k, n* T ,<ri 


where 


<j> • k r i = cb i \ ' = I <{>< k . r \ i d\ = \ 1 — *>\p < - ^ i 


«t>^k, r) = y ) = 2 j \<{><k F r, \ ) d\ - 2 y ~ “ l — 1 1 -*- y i exp i — \ ) I 
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Equation (A. 14) may be written more compactly :n the form 
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Analysis of SAR spectra 


In this section we carry out an mtercomparison of SEASAT SAR image spectra with 
mode! hmdcast wave spectra and SAR image spectra computed from the hmdcast wave 
spectra Most of the analysis of this section is based on a selection of 32 digi tally analyzed 
SAR scenes, listed in Section 5. 1 . However, we present in Section 5.2 also a brief review of a 
previous analysis of a significantly larger set of optically processed SAR images for the same 
SEASAT orbits (Lehner, 1984) 

Theoretical SAR spectra were computed using the linear transformation 
approximation described in Section 4 2 (Section 5 3) and by Monte Carlo simulations using 
the full nonlinear mapping relations described in Section 4 3 (Section 5 5) In addition, 
emoirka! transfer functions were determined by fitting the SEASAT SAR image spectra to 
the hindcast model spectra using a power-law frequency dependent transfer function for 
each wave propagation direction (Section 5 4) The conclusions drawn from these various 
analysis techniques are summarized m Section 5 b 

5 1 SAR scenes analyzed 

A total of 32 digitally processed 64km x 64km SAR scenes from 7 orbits were 
spectral analyzed To test for a possible dependence on the SAR incidence angle, 22 of 
these scenes were selected as pairs from the near and far edges of the 1 00 km wide 
SAR swath The variance spectrum was computed at 25 m resolution (50 m Nyquist 
wavelength) to the full image wavelength limit and then subsequently averaged m the 
wavenumber domain with a top-hat weighting function extending over 5x5 
wavenumber components (yielding 50 degrees of freedom) 

The geographical positions and relevant characteristics of the scenes are given m 
Fig 5 1 and Table 7, respectively Oroits 757, 762, 785, and 791 correspond to the same 
meteorological period, while the remaining orbits correspond to essentially 
uncorretated meteorological situations It is apparent from the nonlinearity ( velocity 
Punching ) parameters and principal wave propagation directions listed m the table 
(inferred from 'he model hmdcast spectra) that the scenes cover a wide variety of 
imaging situations 
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2 Gener al comparison with SEASAT SAR image spectra 

A detailed comparison of the global wave hmdcast with the 32 digitally processed 
SEASAT SAR images, including a discussion of SAR ocean wave imaging problems and 
(.he nonlinear surface wave-SAR image spectral mapping relations, is given in the 
following sections 

For a first overview we show in this sec :on, in Figs 5 2 - 5 5. a comparison of the 
principal wave propagation directions inferred from the variance spectra of a larger 
number of analog processed SAR images (Lehner, 1984) with the Custer diagrams for 
the total sea, windsea and swell of the wave model hindcast. The SAR wave 
propagation rays (panels a) represent the propagation directions of the peak 
wavenumbers of the SAR spectra. The principal wave propagation directions inferred 
for the SAR image spectra show some general agreement with the wave model 
hindcast However, a number of systematic deviations, to be discussed later, are also 
dearly evident: the SAR tends to detect long swell components more readily than the 
short windsea components, the SAR wave propagation direction is often strongly 
rotated towards the range direction (the SAR look direction orthogonal to the satellite 
flight direction), and range travelling waves are more easily detected than waves 
travelling in the azimuthal (satellite flight) direction 

One of the difficulties to be faced in using SAR wave data m a combined wind and 
wave data assimilation system will dearly be to decide whether discrepancies which 
are found between the SAR image spectra and the model wave spectra are due to the 
SAR imaging process or to errors in the input wind field (or the wave model itself) In 
-* the situation is considerably simpler, as alttmeW;' wave 

height measurements are relatively free of distortion and calibration uncertainties 
The effective use of SAR wave data will therefore be critically dependent on the 
development of reliable mapping relations defining the transformation of the surface 
wave spectrum to the SAR image spectrum. The inverse mapping, i e the transformat- 
ion from tne SAR spectrum to the wave spectrum, has as yet not been addressed The 
inverse transformation may be exoected to be non-unique and orobably singular 
Inversion techniques will need to be developed, using general inverse modelling 
methods, employing iterative corrections of the forward mapping relation. These 
questions will be discussed m more detail in Section 6 it .s shown there that the SAR 
inversion problem may not need to be solved explicitly but can be imbedded m the 
more general inversion problem of determining the wind field modifications required 
to correct errors between predicted and observed wave data m a combined wind/wave 


32 


data assimilation system. However, it can be concluded already from Figs 5 2 - 5 5 that 
SAR image spectra clearly cannot be interpreted directly but only by comparison with 
5AR spectra computed from a wave model forecast, and that the assimilation of SAR 
spectra in models will require the development of rather sophisticated inverse 
modelling techniques 

5 3 Comparison with linear theory 

Linear theory should be applicable if the nonlinearity parameter C = C’ nax ro.s-H 
^ l.(C is proportional to the ratio of the rms azimuthal displacement induced by the 
long wave orbital velocity to the dominant wavelength, cf. Section 4 2). Table 7 lists a 
number of cases for which this criterion is satisfied. 

However, m practice it was found that linear theory gave unrealistic SAR spectra m 
all cases, analyzed, including cases with quite small values of C Figs 5 6. 5 7 show two 
examples, corresponding to a strongly nonlinear situation (scene 15601, C = 3 6) and 
the most linear 1 case of the set (scene 11004.C = 0.6) The image spectra are plotted 
in polar coordinates as frequency direction spectra F f.Q i, in accordance with the usual 
wave spectrum format, rather than as wavenumber spectra E<k. = F'f. 9; f!2k‘ 

In both cases the theoretical SAR spectrum shows a generally monotonic increase 
of the spectrum with increasing frequency, without a distinct spectral peak as found m 
the observed SAR spectra The origin of this behaviour is readily identified The 
theoretical velocity bunching SAR MTF is proportional to fl (cf. Section 4 2) Thus 
F.SARifc&HFtmvef’Q) ~ t* Smce the wave model hmdcast spectra typically exhibit a high 
frequency dependence - f-*. the SAR spectra increase monotomcally as f 2 for high 
frequencies. A SAR peak is visible at low frequencies only if the wave spectral peak is 
sufficiently sharp that it overrides the f* factor, i e. if 


. tvave 

/ — — 


i/', 9) < -F,F if' 9) 

unv 


(5 1) 


for some frequency on the rear slope of the peak This is seldom observed, even for 
swell spectra (Bruning et al , 1988. computed such cases for theoretical swell SQectra 
represented by the JONSWAP formula with very high peak enhancement factors 
V = 10 ). 
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FiS 5.2 Ray diagram (panel (a)) constructed by lehner(1984) from peak wave- 

numbers of 5EASAT SAR image spectra along orbit 757, August 18, 1978, 

22:40, with superimposed Custer diagrams of wave model hmdcast for total 
wave field (panel (b)), windsea (panel (c)), and swell (panel (d)) 
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Ray diagram (panel (a)) constructed by Lehner (1984) from peak wave- 
numbers of SEASAT SAR image SDectra along orbit 762. August 19. 1978. 

06 40. with superimposed Custer diagrams of wave model hindcast for total 
wave field (panel (bj), wmdsea (panel (0). and swell (panel (d)) 












F| q 5 4 Ray diagram (panel (a)) constructed by Lehner ( 1 984) from peak wave- 

numbers of SEASAT 5AR image spectra along orbit 785. August 20. 1978, 

21 40. with superimposed Custer diagrams of wave model hmdeast for total 
wave field (panel (b)). wmdsea (panel (c)), and swell (panel (d)) 
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Ray diagram (panel (a)) constructed by lehner (1984) from peak wave- 
numbers of SEASAT SAR image spectra along orbit 1 359, September 30, 1 978, 
01:17, with superimposed Custer diagrams of wave model hindcast for total 
wave field (panel (b)), wmdsea (panel (c)), and swell (panel (d)) 












S^CCTftUM SAJI SPCCmuu ORBIT 783 SCENE MO 04 TRAMSFO8MC0 UOOCt SPECTRUM 



Observed SAR spectrum. wave hmdcjvt spectrum and theoret 
imaged SAR spectrum (or SAR scene I 1004 (linear. #=06! 



The failure of linear theory to reproduce the observed shape of the spectrum is not 
m contradict on with the theoretical range of validity of the theory Linear theory- 
breaks daw- for short wave lengths of the same scale as the rms azimuthal dispiace- 
ment The azimuthal displacements associated with these components become 
smeared out by the superimposed displacements of the dominant long waves, which 
shift the small scale velocity bunching pattern back and forth relative to the phase of 
the short waves The peak occurs m the SAR spectrum in the region of transition from 
the linear to the nonlinear regime On the low frequency side of the peak, linear 
theory is still applicable. However, to test the linear theory quantitatively in this 
regime would require an investigation with calibrated images, which was not 
attempted in the present study (This would not require an absolute calibration of the 
SmR , as the SAR image is defined as the modulation image with respect to the 
background white noise clutter image. However, it would require a more detailed 
investigation of the reference clutter spectrum than was undertaker? here,) 

5 4 Empirical transfer function fits 

It is of interest to investigate whether the SAR image spectrum can be related to 
the surface wave spectrum through an empirical linear transfer function Although we 
have seen that the purely linear theory is generally inapplicable for the full spectrum, 
it has been argued that the principal effect of the orbital motion is the velocity spread 
term, whicn can be represented by a linear az>muthal cut-off filter function (cf Beale, 
1983) Also, the exact nonlinear transformation expression (4 2) derived m Section 4 
has a quasf-Imear (although non-diagonal) form 

Accordingly, empirical MTFs were determined for a set of 26 SAR soectra (six of the 
cases listed in Table 7 could not be analysed as the scenes were received too late) by 
least square fitting the free parameters A and n m a power-law ansatz of the form 

f <:») = w)»i - r”* 

v\K ' j.' i 

where the reference frequency f 0 was held fixed at 0 i Hz 

Figs 5 8, 5 9 show two typical examples of the observed and fitted SAR spectra, 
together with wave model hmdeast spectra tnd the directional dependence of the 
fitted parameters AiH>, mti/ In nearly all cases quite satisfactory fits were obtained 
The exponents ni 0; of the power laws were negative throughout, typically varying 


F ~ F ' r , h *■ ISO i (5 2) 

cm* v*’t» 
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f between - 5 and - 1 0, in contrast to the predictions of the purely linear theory This is 

presumably due to the smearing of the modulation of the shorter waves by the 
supem.iposed azimuthal displacements induced by the longer waves, as discussed 
above, which produces an attenuation of the shorter waves This velocity bunching 
effect generally produces a significant attenuation of the spectrum already in the 
intermediate scale range shortly to the right of the wave spectrum peak, whereas the 
azimuthal cut-off due to the velocity spread term becomes effective only at higher 
wavenumbers (cf. Section 4). 

Pig 5 1 0 shows the dependence on wave propagation direction 0 relative to the 
satellite flight direction of the mean values and the standard dewri on of the 
exponent n($) for the 26 cases studied The ratio of the standard deviation to the mean 
(the scatter index) is clearly too great for these empirical MTF fits to be applied in 
oractice as a method for predicting the SAR spectrum. However, they do summarize m 
a single form the mam characteristics of the shape of the SAR spectrum m the energy 
containing range relative to the wave spectrum. 

A closer inspection of the data indicated that the scatter seen m Fig. 5. .0 cannot be 
significantly reduced by a stratification of the data with respect to the nonlinearity 
parameter Cor other spectral parameters We conclude from this and the previous 
section that a reliable interpretation of SAR image spectra must be based on realistic 
com putations of the full nonlinear mapping of the two-dimensional wave spectrum 
into the SAR image spectrum 

5.5 Monte Carlo computations 

The numerical computation of the closed nonlinear integral transform expression 
(4 31) relating the surface wave spectrum to SAR image spectrum could not be 
implemented within the present project. This is planned for a follow-up study project 
it is anticipated that the technique, when implemented, will be more accurate and 
efficient than the Monte Carlo methods used hitherto For the present project 
numerical computations of the nonlinear spectral transformation based on the same 
general SAR imaging expression (4.20) were carried out using the existing Monte Carlo 
technique (Alp^rs and 8runmg, 1986, Brunmgetal , 1988) 
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spectra, together with wave spectrum 









Observed and empirically fitted SAR spectra, together with fitting parameters 
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Fig. 5.11a-| shows a comparison of observed and computed SAR spectra, together 
with the wave modei hindcast spectra used for the Monte Carlo computations. The 
selected set or SAR spectra covers a representative cross-section of the 26 cases studied 
(cf. Tabel 7). The examples include range and azimuthally travelling waves, weakly and 
strongly nonlinear imaging regimes and single and multiple peaked spectra. 

The agreement between the observed and simulated spectra is generally quite 
remarkable. The significant distortions of the original wave spectrum appearing tn the 
SAR spectrum, m particular the pronounced azimuthal cut-off, is well reproduced in 
the simulations. The propagation directions ana wavelengths of the peaks in the 
simulated and observed SAR spectra also agree quite well A somewhat earlier 
azimuthal cut-off is generally found in the simulated spectra than in the observed SAR 
spectra This could indicate a rather too high energy level in the high frequency region 
of the wave model spectrum. 

The sequence Fig. 5. 1 1 , panels a, b, c. and d, give a good example of the change in 
character of the spectral distortion as one progresses through a region of azimuthally 
travelling waves to a region of range travelling waves along orbit 757 (cf Figs 5 1 and 
5 2). The discrepancies seen in Fig. 5 2 between the mean propagation directions of 
the model wave hindcast and the propagation rays inferred by Lehner (1984) from the 
peak wavenumbers of the SAR spectra are clearly due to the azimuthal cut-off in the 
SAR spectrum of the predominantly azimuthally travelling waves along the north- 
eastern segment of the orbit 

Although it is encouraging th< t the nonlinear SAR imaging theory is able to 
reproduce the principal features of the observed SAR spect r e. Fig 5 11 aiso illustrates 
the problems one faces in assimilating SAR image spectra m wave models: the loss of 
information for azimuthally travelling waves, the 180° directional ambigu •ty (Fig 5 11, 
panel j, gives a good example of two wave fields propagating in opposite directions 
which are merged together in the SAR spectrum), and the pronounced distortion of 
the spectrum None the less, the potential information content of the two-dimensional 
SAR spectrum is clearly greatly superior to a single scalar number, such as the 
significant wave height, for example, or even a one-dimensional frequency spectrum, 
if effective techniques for the appropriate exploitation of this information can be 
developed (It should be noted in this context thar SAR wave image spectra are inde- 
pendent of absolute SAR calibration uncertainties: The spectra are 'self-normalized' by 
the clutter background level.) 
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Examples of observed SAR spectra, wave model hmdcast spectra and SAR 
spectra derived from the wave spectra by Monte Carlo computations 



SAH SPtCIRUM: OIKIIT 7 «2. SCfNT 323. SIMULATED SAR SPECTRUM 
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Examples of observed SAR spectra, wave model hmdcast spectra and SAR 
spectra derived from the wave spectra by Monte Carlo computations 
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Examples of observed SAR spe tra. wave model hmdcast spectra and SAR 
spectra derived from the wave spectra by Monte Carlo computations 
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model hmdcast spectra and SAR 
' Monte Carlo computations 





5 6 Conclusions 


The principal conclusion to be drawn from our SAR studies may be summarized as 

follows: 

( 1 ) The linear SAR imaging theory is not able to reproduce the shape of the observed 
SAR image spectrum, even for low energy sea states. This is due to the unrealistic 
amplification of short waves by the linear theory, which results in a monotom- 
cal ly increasing SAR spectrum In realuy, the higher wavenumbers (often 
beginning shortly after the spectral peak) are attenuated by nonlinear inter- 
action* with the motion induced azimuthal displacements of the longer waves 
For low energy sea states with small nonlinearity parameter C n ' ax , the observed 
peak in the SAR spectrum represents the transition from the linear to the 
nonlinear imaging regime 

(2) In nearly all cases, the observed SAR spectra can be well reproduced by applying 
individually fitted, directional dependent power law MTFs to the hmdeast wave 
spectrum However, the parameters derived from individual fits show 
considerable scatter from case to case Thus this method does not appear suitable 
as a general empirical tool for relating wave spectra to SAR image spectra 

(3) SAR spectra computed from hmdeast wave spectra using Monte Carlo simulation 
techniques showed a remarkable agreement with observed SAR spectra Smce 
the nonlinearly distorted SAR spectrum can differ significantly from the original 
wave spectra, and the distortions depend critically on the detailed structure of 
the wave spectrum, a meaningful interpretation of SAR wave image spectra can 
be based only on numerical computations of the exact nonlinear transformation 
relations 

(4) Starting from the well tested theory describing the nonlinear mapping of the 
ocean wave scene into the SAR image plane wh.cn was aopiiec m the Monte 
Carlo simulations, a closed nonlinear integral exDression was derived for the 
transformation of the wave spectrum into the SAR image soectrum Since Monte 
Carlo computations are too time consuming to be aooued operationally tor the 
analysis of ail ERS- 1 SAR wave mode data (approximately I minute comoutmg 
time per spectrum on a CRAY X-MP), the numerical implementation of the new 
expression should be pursued and tested with high priority as a possible alter- 
native for application to ERS- 1 


(4) The oroblem of the inverse transformation from the SAR image spectrum to the 
surface wave spectrum has not yet been addressed Since different back scattering 
elements can be mapped on to the same image pixel, the inversion is in general 
non-unique and singular The inversion problem must therefore presumably be 
treated locally in the form of a perturbation about a given forward 
transformation. This implies that the operation of a wave model, which can 
provide the first guess wave field needed as input for the forward 
transformation, is an essential pre-requisite for the extraction of useful wave 
information from SAR image spectra The optimal implementation of the 
inversion procedure will depend on the form m which the SAR wave data s 
assimilated m the overall wind and wave data assimilation system in development 
for ERS-1, which is discussed further in the following section 
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6. Data assimilation 


As outlined in the introduction, the present study represents a first step towards the 
planned implementation of a comprehensive wind and wave data assimilation system for 
ERS-1 A characteristic feature of this system is the joint analysis of both wind and wave 
data, with consideration of the cross coupling between wind and wave information in the 
different ERS-1 sensor systems and in the assimilation procedure itself In this study only a 
small sub-set of the questions to be addressed in the implementation of the full system will 
be considered. In the following sub-section 6 1 we d.scuss briefly the general problem of 
assimilating wave data in wave models and the information which can be extracted from 
wave data, via a wave model, for the wind field, in Section 6.2 we present a Simple example 
of the improvement in a wave forecast obtained by assimilating SEASAT altimeter wave 
height data 

6.1 G eneral considerations 

In contrast to atmospheric data assimilation, which has a fairly long history 
beginning with the introduction of numerical weather prediction, wave modellers 
have only recently begun to address the problem of assimilating wave data m wave 
models (cf studies by Thomas (1988), Komen(1985), Janssen etal (1988) This is partly 
due to the lack of sufficient wave data in the past to significantly affect the wave field 
analysis Here ocean satellites will, of course, radically change the situation But wave 
modellers have also not viewed data assimilation with the same urgency as atmo- 
spheric modellers, since wave forecasting, as opposed to weather forecasting, is not 
essentially an initial value problem, but a boundary value problem, the wave field is 
forced by the giver, wind field, and the initial wave field at some far past time when 
the model was first spun up is irrelevant. Wave data is used not to initialize the wave 
field, but simply as a running correction or update of the model simulation. In fact, it 
could be argued that the mam interest m wave data assimilation is not so much in the 
improvement o* the wave forecast but m the opportunity it provides for improving the 
wind field analysis. For the wave field is a highly sensitive indicator of the Quality of 
the wind field, as demonstrated m Sections 2 and 3. 


T he difference m the nature of wave data assimilation compared with standard 
meteorological data assimilation techniques carries a number of implications (cf 
Hasselmann, 1985): 
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(i) The continuous assimilation of wave data during the running of the model creates 
no dynamical shock' problems. The local update information simply propagates 
into the forecast in the same way as a change in the external wind forcing 

(li) If the modified wave field consists primarily of swell (as m most regions of the 
ocean, cf Section 2), the updated swell will continue to propagate at its corrected 
value throughout the rest of the forecast: the correction sticks' However, if the 
wave field represents a wmdsea, and the (incorrect) wind field responsible for the 
wave field error is not modified simultaneously with the wave spectrum, the 
updated wave field will rapidly relax back to its original incorrect state in response 
to the uncorrected wind field. 

(iii) Wave data is typically incomplete. The model requires the field of two-dimensional 
wave spectra, which are estimated today at only relatively few locations with 
directional buoys ERS-1 will provide data on a global scale, but only in the form of 
the altimeter significant wave heights and azimuthally clipped, nonlmearly 
distorted and 180° ambiguous two-dimensional SAR image spectra. The 
assimilation procedure must therefore address the problem of translating this 
incomplete information into optimal spectral updates. 

(iv) The region of influence of a measurement at a given point cannot be taken as 
approximately constant, as in the atmospheric case, but depends strongly on the 
type of spectrum and the spectral component. Swell which was generated in a far 
distant source may be associated with a region of influence of several thousand 
kilometers, while wmdsea generated by a local intense storm will have a much- 
smaller characteristic correlation scale of the dimension of the wind field We show 
in the following that the proper treatment of this problem automatically follows 
from a consistent coupled wind and wave data assimilation approach. 

6 2 Optimal wave data assimilation 

The general task of wave data assimilation i$ to optimally modify the total wave 
field F f t\ 0. x. t, predicted by the model up to the timer for all points x on the global 
ocean (or some selected area, in the case of regional applications) on the basis of a 
finite set of measurements ./, =<p,iF . i - / n 


The data d. can represent single numbers (e g. H<) or a full SAR spectrum, for 
example, it will be assumed that the wave data are assimilated continuously, so that 
we assimilate only wave data given at the time t. 

Optimal' implies here that the modified wave field should be consistent both with 
the data d, and thedynamica 1 constraints of the wave model In addition, the wave 
field is normally required to satisfy some additional aesthetic' quality or plausibility 
criterion, such as smoothness or minimal least square deviation from the predicted 
wave field. The dynamical model constraint implies that a consiste.it modification of 
the wave field can be achieved only if the wind field (at some earlier period) is also 
m xdified The quality criterion will therefore be extended to include the requirement 
th st the modification of the wind field should also be kept as small as possible 

6.2 1 General formalism 

Satellites and conventional wave instruments yield various types of wave data 
djx p , t p ) at a given set of measurement positions x., and times t p which are related 
to the two-dimensional wave spectrum Fix p . t p . through functions 

d = DtF) < 6 ’> 

The goal of wave data assimilation is to modify the wave spectrum F< x n . t p/ 
predicted by the model m such a manner that the wave data error function 

e = ^ >d - d ) M id - d ) ., 

;t — •■'.’)! (6 2 ) 

defined with respect to some suitably chosen (normally diagonal) positive definite 
matrix is minimized, while maintaining the internal dynamical consistency of 
the overall wave field. 

The latter condition requires that the wave spectrum modifications AF should 
be induced through modifications An. in the driving surface stress field, 

AFlx . M = j P dt'\ dx'Cix. v . / . x'.f'Uu.ix’.n (6 3) 

• - c J J - } 

or A F = G< Au*; 




Here G( x pj t p ; x' t*) represents the (Green function) response of the wave field 
to small perturbations of the surface stress field It is assumed in (6 3) that the wave 
field perturbations are small, so that the relation between IF and Au, can be 
linearized. 

In general, the modified stress field will not be uniquely determined by the 
condition e, L , = minimum. We therefore add an additional constraint: we seek a 
modification of the wind stress field which minimizes n :u , while at the same time 
keeping the changes m u. as small as possible This can be achieved by minimizing 
the net error function (taking now for simplicity diagonal weighting functions w t 
and w n ) 


^ iv ( (Id + d - d) ) j dt' dx'\ Au # (x\ H ] -w f*) 


with respect to Aiu, where A d r = D, 

The minimization of e is to be carried out continuously as a function of t, 
equation (6.4) containing at any time t all measurement points t p with t p s t In 
practice, the minimization of e at a given time t will be able to benefit from the 
previous minimization at time t - A t, where A t is the time increment for the 
insertion of new data (in the simpler example considered in the following sub- 
section, A t is the 20 minute propagation time step of the global wave model) 


The full problem of minimizing e as given by equ (6 4) is dearly a formidable 
task. First, the Green function G has to be determined This can be avoided by 
applying the adjoint technique, which is based on iterative integrations of the 
wave transport equation and its adjoint equation (cf Thicker and Long, 1988) 
Secondly, e has to be minimized by some iteration procedure with respect to the 
field u ,'x. r depending on two space and one time dimension For a system with a 
lar ge number of degrees of freedom, this will normally require many iterations 
Since the integration of the adjomt wave transport equation alone requires a 
computational effort comparable to the integration of the wave model itself, the 
adjoint technique does not appear very promising without introducing additional 
simplifications. We shall follow instead the direct approach of inverting the wave 
transport equation and develop a very simple approximation for G 


6 2.2 Approximate treatment of the Green function G 


The relation between a perturbation bu, of the surface stress field and the 
resulting perturbation induced in the wave spectrum is obtained by perturbing the 
basic wave model transport equation (cf WAMDIG, 1988) 

(i +vV j F = s ^ s ^ Sj (6 5) 

where v denotes the group velocity and S.,. S^i. $&, represent the input, nonlinear 
transfer and dissipation source function, respectively (For simplicity, we here used 
cartesian rather than the spherical coordinates of the global WAM model.) 

One obtains 


/ d rifi 

( ■+■ vV \fjF — A6 F +■ F — ■ .5u, 

' * du. 


(6 6 ) 


where the linear operator A represents the functional derivative of the nonlinear 
source function .5 . es . t = S t „, * S,„, excluding the wind input term 


and the second term on the right-hand side represents the derivative of the linear 
input source function S. n = pf with respect to u. 


The operator A has been discussed in WAMDIG (1988), since it arises in the 
implicit source function integration scheme used m the WAM model. 8ecause of 
the nonlinearity of the source functions and S ts , .t is in general non-diagonal 
his implies that a n- function perturbation 6u introduced at a given grid point ana 
time generates not omy wave spectrum perturbations propagating away from the 
disturbance point with their appropriate group velocities (as described by the 
diagonal propagation operator on the left-hand side of equ (6 6)), but also a 
secondary perturbation field produced by scattering from the primary rays into a 
continuum of non-focussed secondary rays 


This greatly complicates the computation of the Green function The same 
scattering process arises and is properly computed, of course, in the integration of 
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the fully nonlinear wave model. However, to avoid having to carry out a similar 
computation to determine the Green function, the matrix operator A may be 
simplified by retaining only the diagonal part A . The same approximation was 
used in the implicit integration scheme of the WAM model and was shown to be 
quite satisfactory for that application. In this case, equ. (6 6) can be immediately 
integrated to yield the Green function 


\F = C(AiO = 



?( d JL 

‘ air 


Air 


exp 


[ dfUr'.n 

J t' 


(6 7) 


where the integrations with respect to r'and t" represent path integrations along 
the individual wave group propagation paths 


x' - X ^ V -\t' - t) 


(6 8 ) 


X” = .-I r- f) 

terminating in x and x*, respectively. 

An inspection of (6 7) reveals that for any given spectral component G has a 
pronounced maximum at a position along the ray corresoonding to the last pomt 
at which the wave component received a significant input from the wind, i e m 
the transition region from wmdsea to swell For earlier times, the impact of a 
variation 6u. is lost through the interaction with the rest of the spectrum, 
expressed by the large (wmdsea) value for V In the swell region, a*ter the wmdsea 
has been transformed into a swell component, the exponential damping vanishes, 
but the wind input factor ofl/du w itself also vanishes. The most effective region for 
modifying a spectral component is therefore essentially identical to the 'region of 
origin of the wave component as defined heuristically through the 'wave age i 
(cf 8ooi] and Holthuijsen, 1987) 

: o ; resonable first approximation the Green function may therefore be 
represented as a S-function: 

Gix, t . ; x\ t') = Si x* - x -h vtJou' - t - a A 

where A = .Ai k, x, l ; is a weighting factor determined by the speed of transition 
from windsea to swell 
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The relatively simple structure of the approximate Green function and the 
reduction of its effective 'region of dependence' to a geometrically weii defined, 
limited region governed by the wave age suggests further simplifications in the 
next step of determining the A//, field which minimizes*:. 

6.2.3 Implementation 

We nave established linear relationships between the three sets of data 
consisting of the perturbations of the wind stress field ou. ^ ( ^ a j r perturbations of 

A 

the wave spectra 8F * r^and perturbations<v</ ( -di> = of the deviations 
between predicted and observed wave data 



a 


(6 10 ) 


D 


— \a r a 


(6 11 ) 


where Af aa 1 ' represents the Green function (equ 6 9) and .V/ ‘ represents the 
derivative of the predicted data d t with respect to the predicted spectrum f The 
indices of the vectors 4> a , D v and n a run through the entire data set involved in a 
given data assimilation cycle, including the spatial and temporal dependence and 
the individual components of the fields (e g. the set of spectral components in the 
case of wave or SAR image spectra, or the two vector components m the case of 
the wind stress vector) 


In terms of the present notation, the variational problem (6.4) we wish to solve 
takes the form 


, - V 


w \D 


u - min 

■i 


(6 12 ) 


where 


D 



(6 13) 


represents the deviation between the predicted and observed wave data for the 
unperturbed wave spectra (prior to data assimilation) and w lU w x are weighting 
factors 
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Substituting the linear relations (6 10), (6 1 1) into (6 12) we obtain 


e - ^ w ( d° - ^ .w;; n a ) * w it ^ = 


min 


(6 14) 


.W ,a ’ = V ,v/‘ 21 .Vf' 1 ‘ 

•ti — n oq 


(6 15) 


The minimum of (6. 12) is given by the solution of the linear system of equations 


^ a/' 4 ; 

<l|l 

U 



(6 16) 


where 


•V/' 41 = ^ w .Vf' 31 

OP — — .Q VP 


w o , 

I ) up 


(6 17) 


and 


C = ^ w D u M' 3 ' 

u — * v v vn 

in practice, the minimum of £ may be obtained more rapidly by applying direct 
minimization methods (e g conjugate gradient techniques, cf. Navon and Legler, 
1987) to the relatively simple quadratic form (6 12) rather than solving the system 
(6 14) 


The solution of the minimization problem yields the entire wnd stress 
correction field 6u,, but the corrected wave field only at points where measure- 
ments exist. (The matrix multiplication (6 1 5) runs over only those indices for which 
1 e on| y those spectral values which impact the data values <?,-.) To up- 
date the entire correction field for the wave spectrum the Green function relat.on 
(6 10) must be applied for the wave spectra at all grid points. 

We note that the aoproach outlined here avoids the need for an explicit 
inversion of the wave spectrum — SAR spectrum transformation relation The 
inversion for SAR wave data is considered automatically within the framework of 
the complete inversion problem. We note furthermore that the derivative matrix 
^va ~‘ for SAR im age spectral data was derived in Section 5 in the course of the 
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derivation of the general nonlinear spectral transformation relation for SAR image 
spectra 

Since effective numerical procedures have been developed for minimizing 
quadratic forms of the type (6. 1 4) for systems with several thousand degrees of 
freedom (cf Navon and Legler, 1986), it appears feasible to implement the 
approach outlined here in operational framework. 
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As a first step towards the development of a complete wind and wave data 
assimilation system we have investigated the simpler decoupled problem of assimilating 
wave data into the wave mode! without modifying the wind field This differs from the 
approach of Thomas (1988) and Janssen et al (1988), who modify both the wind and wave 
field. As outlined above, our ultimate goal is to modify only the wind field directly, the 
wave field changes being tied to the wind field modifications through the wave model The 
present decoupled study was earned out as a reference to establish the direct impact of the 
wave data on to the wave model. It may also be interpreted as an implicit wind field modi- 
fication exercise in which only the past wind field (which has no influence on the current 
wove field development) is altered. The study provides insight into the characteristic region 
of influence of wave measurements, and the typical relaxation time of new information 
inserted into the model The assimilated data were the global SEASAT altimeter wave 
height data obtained during the month of August, 1978 This month was chosen as it 
contained relatively few data gaps 

The data were inserted into the wave model at each propagation time step, t e every 
forty minutes This corresponds to an orbit segment of about 1 7,000 km, or up to 85 data 
points separated by 200 km. For each wave height data point the model wave spectrum was 
modified over a finite set of surrounding model grid points Experiments were carried out 
with different scales L for the region of influence 

Within a region of influence, an updated wave spectrum F nt ?, v 'f. O' was formed by 
applying a correction factor c to the predicted spectrum F pri >j f. 0,, 

F i/, 0) = cF »/. 8) ^ ^ 

neic pr*»d 


where 


c - 


(i. 

lA<J>r 

iA.\r 

’ 2 


L 2 

\ 


l for lAif) 1 > <>r QAi > 


> 2 


) ) for lA»{>l < L . IA\ < L 

' ' p \ 


and 
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Here H s alt . H s mo ^ denote the altimeter wave height and the interpolated model 
wave height at the measurement point, respectively, and A4>, A\ represent the distances »n 
degrees latitude and longitude respectively, of the model grid point from the measurement 
point (cf. Fig. 7 1) Close to the measurement point, c * cmax* and the model spectrum is 
adjusted to reproduce the observed wave height. The correction factor relaxes bi-lmearly to 
unity at the edge of t ne rectangular region of influence 

Assimilation experiments were carried out with - 9. 9 and 15, 15 degrees, 

corresponding to 3 x 3 and 5 x 5 grid points in a region of influence The best results, 
shown m Figs 7 2 - 7 5, were obtained for the case L = 15 This scale is compatible with the 
time-space sampling density of the satellite and ensures that essentially all points on the 
model grid have received significant impact from the altimeter wave height measurements 
within a 3 day sampling penod 

To remove spm-up effects from the August analysis pertod, the assimilation run was 
started on July 28, the assimilation run itself being restarted from the run without 
assimilation beginning July 7, 1978 The GLA surface stress field was used as input 

Figs 7 2 a - c show the mean wave heights for August, obtained with the assimilation 
run A comparison with Figs 2 2 o - 7 4 b and 2 5b showing the previous results for the run 
without data assimilation and the sigm ficant wave heights derived from the S6ASA^ 
altimeter demonstrates the significant improvement achieved through data assimilation 

The relaxation time of the assimilation proceoure was tested m a further experiment in 
which the assimilation process was switched off in the middle of the (3-pomt) assimilation 
run, on August 1 5 The characteristic relaxation time for the global model bias error (the 
global mean deviation of the model and altimeter wave height) is of the order of 6 days 
(Fig. 7 3b) in this oeriod the error nas grown oack to the same magnitude as in the original 
hmdcast without data assimilation 

We note that the characteristic relaxation times found for the oresent ass.miiation 
switch-off experiment are generally larger than those found by Janssen et ai (1988) m a 
similar switch-off experiment. This isdue to the larger region of influence chosen in the 
present assimilation scheme. If the region of influence is smaller than the orbit spacing 
corresponding to the satellite repeat period (or characteristic time scale of weather 
patterns, »n the case of wmdsea) only part of the model grid is updated during the 


assimilation procedure The propagation of the non-corrected wave field from the missed 
regions into the neighbouring, corrected wave field regions then degrades the update m 
the con ected regions on the relatively short time scale of the wave propagation time from 
non-updated to updated regions 
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7 2c Mean significant wave heights for model hmdcast with assirr Hated altimeter 
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8. Conclusions 


Our investigations of SEASAT altimeter and SAR wave data using a global wave model 
driven by wind fields derived from assimilated scatterometer data have provided valuable 
insight into the quality and application of ocean satellite wind and wave data. This SEAjAT 
exercise is particularly relevant for ERS-1 , which will be carrying a similar suite of sensors 
We may summarize the principal conclusions of this study, which we regard as a first step 
towards the development and implementation of a comprehensive wind and wave data 
assimilation system for ERS-1, as follows; 

1' Wave models represent an excellent data validation tool, providing a sensitive natural 
consistency check of wind and wave data The shortcomings of the GlA 6 hourly 
averaged surface stress fields derived by data assimilation from scatterometer and 
conventional meteorological data were immediarply rpupalpn n\j comparing rhp 
hmdcast wave heights with altimeter wave heights This technique is particular 
valuable in data sparse regions such as the Southern Ocean, when the altimeter wave 
height data may represent the only independent verification of the wind field analysis 

2) A number of further data validation tests and wave hindcast intercomparisons with 
other wind field products suggested that these shortcomings were not directly related 
to the wind field data assimilation method as such, but resulted from the time 
averaging involved m producing the stress fields But it »s also poss.b^e that 
inconsistencies m the scatterometer algorithm (Woiceshynetai . 1987. Anderson et al . 
1987) contributed to the discrepancies A combination of the data validation 
procedures described in Anderson et al (1987), Janssen et al (1988) with those in the 
present report would provide a valuable operational system for real time monitoring 
of the performance of the wind and wave sensors of ERS-1 

3) T he direct assimilation of altimeter wave height data in wave models has a significant 

positive impact on the skill of the wave forecast This applies particularly for the swell 
field The average me of : t ^y T v^^p jjnc^atp i n fp r\ f f ih a. o. r dp r n f c, * 

days 

4) For winosea:, wave data assimilation has little impact on the quality of the longe" term 
wove forecast unless the wind fieid is moamea togetner with the wave held By 
improving the wind field, this requirement also increases the skill of the weather 
forecast and provides improved air-sea flux fields required for climate research 
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5) The development of a fully interactive wind and wave data assimilation system 
represents a formidable task, since wave field errors imply errors in the wind field not 
only at the same time level, but also at early levels which can extend back one or two 
weeks. A general variational approach for constructing simultaneously optimized 
wind and wave fields was developed The models and numerical techniques for 
implementing such a system are now largely available, and it is recommended that 
work in this area is pursued with high priority in order to provide a viable operational 
system in time for the launch of ERS-1. The modifications of the wind field induced by 
wave data in such a system would represent a quantitative, constructive version of the 
wind and wave data mutual consistency tests summarized in 1) above 

6) SAR wave image spectra normally represent strongly nonlmearly distorted trans- 
formations of the wave spectra. This applies for both wmdsea and swell spectra SAR 
image spectra cannot be usefully interpreted unless supported bv commutations of t h e 
nonlinear transformation from the wave spectrum to the SAR spectrum . This implies, 
m practice, that SAR wave spectra can be usea operationally only m conjunction with a 
realistic (third generation) wave model which can provide a good first guess wave 
spectrum This will be available only man operational forecasting environment, i.e. in 
weather forecasting centres. 

7) Monte Carlo computations of the nonlinear transformation of the wave spectrum mto 
the SAR image spectrum, using the hindcast wave model spectrum as input, were in 
excellent agreement with the 5EASAT SAR spectra for a set of 32 different scenes This 
indicates that the long debated theory of the SAR imaging of a moving ocean wave 
surface is now well understood and can be used confidently 

8) A new closed expression relating the SAR image spectrum to the surface wave 

spectrum was derived This opens the possibility of computing the nonlinear spectral 
transformation directly, without reverting to the pixel-to-pixe 1 mapping of individual 
scene realizations using Monte Carlo simulation techniques The Monte Carlo 
technique suffers from stat^sthC^I^SHa^mjD^ing errors and is r ‘ r * ,|- stvpijV' 

time I” 1 minute CRAY X-MP per spectrum) It cannot be used for the routine ptocessmg 
of the many SAR image spectra generated globally by the ERS- 1 SAR m the wave mode 
Computations of the nonlinear spectral transformation using the newly derived 
relation were not foreseen and could not be included within the framework of this 
study contract, but should be carried out with high priority in a follow-up study to 
determine whether thisappioach has operational potential for ERS- 1 
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9) Methods need to be developed for inverting the transformation wave spectrum $AR 
spectrum. This question should be approached within the general context of imple- 
menting a joint wind and wave data assimilation system as summarized under 5) 
above. The goal of such an interactive assimilation system is to modify only the wind 
field directly, the wind field modifications being chosen such that both the wind field 
and the resultant wave field are optimized with respect to the wind and wave 
observations In this context a SAR image spectrum can be treated in the same manner 
as any other derived wave information, such as altimeter wave heights or one- 
dimensional wave spectra. It follows that the inversion of the transformation wave 
spectrum SAR image spectrum need not be carried out explicitly, the general 
inversion procedure being based on iterative corrections including only forward 
transformations. 

In summary, we have found encouraging results with respect to the potential of wave 
models for the routine validation of wind and wave data and thp a«imilatinn nf motor 
wave data, 'ost of these could be readily implemented in an operational system However, 
we have also identified (and in some cases made some progress m theoretically resolving) a 
number of fundamental issues with respect to the simultaneous, interactive assimilation of 
wind and wave data and t h e application of SAR wave image data These questions will need 
to be addressed with high priority within the next two years if ERS- 1 is to realize the high 
expectation placed by the scientific community in this satellite 
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GLOBAL DATA ASSIMILATION PROGRAMME FOR AIR-SEA FLUXES 

J$G / CCCO Working Group on Air-Sea Fluxes 


Summary 


The JSOCCCO Working G r oup on Air-Sea Fluxes was formed by the Joint 
Scientific Committee (JSC) of the World Climate Research Programme (WCRP) 
and the Committee on Climatic Changes and the Ocean (CCCO) to develop a 
programme for producing, on a continuous basis, global fields of air-sea fluxes 
.making use, in particular, of the new data which will become available with the 
next generation of ocean satellites These fields are an essential input for many 
climate studies, including the major WCRP/CCCO projects TOGA (Tropical 
Ocean/Global Atmosphere) and WOCE (World Ocean Circulation Experiment) 

The goal of the proposed Global Data Assimilation Programme for Air-Sea Fluxes 
(GDAP) is to provide continuous, 6 hourly, 1° resolution, global fields of the air- 
sea fluxes of momentum, sensible and latent heat, water, and sola? ; md infra-red 
radiation beginning 1991 


To achieve this goal, comprehensive data assimilation systems need to be 
implemented at operational global weather centres capable of processing the 
many different forms of input data which enter into to the flux field analysis. 
These include data from oceanographic and meteorological satellites, ships, 
buoys and the weather station network it is proposed that at least one such data 
assimilation system be established in Europe (at the European Centre for 
Medium Range Weather Forecasts) and in the U S (at the National 
Meteorological Center). Additional assimilation systems at other weather centres 
could usefully augment the programme 


The 1991 target date set for the operational implementation of the system at 
ECMWF and NMC corresponds to the estimated start of the operational phase o* 
ERS-1, the first of a senes of planned oceanographic satellites Before tne 
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phase the system wiil be run in a research and validation mode, with 


no direct interaction with operational analyse? and forecasts 
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A number of projects are defined to develop and implement the system within 
this time frame. These include: 


the incorporation of ocean satellite sensor algorithms within the 
assimilation system, 

the quasi real-time provision of calibrated physical sensor data from ocean 
satellites, 

the extension of existing data quality control systems to include the new 
kinds of data available from oceanographic satellites, 
the inclusion of surface wave data and global wave models m the sensor 
algorithms, data assimilation system and flux parameterisation, and 
t in improvement of assimilation procedures for the insertion of two- 
dimensional surface data into three-dimensional fields 

In addition to the near real-time data assimilation programme carried out within 
the framework of the routine forecasting operations, provision should be made 
through an effective archiving and retrieval system for a possible reanalysis of 
the entire input data set at a later time. Although costly, this could prove to be 
desirable after a few years operational experience in order to provide a 
temporally homogeneous flux data set based on the latest improved algorithms 
and assim'lat'on systems 

Proposals are also made for upgrading the in situ ocean oata collection system In 
situ data are important both as valuable input data and for calibration purposes 
Organizational measures are suggested which could Significantly enhance the 
quality and quantity of data from these sources 
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Goals 


The fluxes of momentum, heat, rad.ation and water at the air-sea interface 
drive the ocean circulation and determine the coupling between the atmo- 
sphere and the ocean. A reliable determination of these flux fields is there- 
fore essential for understanding the role of the oceans in climate The goal 
of the Global Data Assimilation Programme (GDAP) is to develop practical 
oceanographic satellite data retrieval and appropriate data assimilation 
schemes which could be applied routinely in producing optimal global air- 
sea flux data sets. 

This activity has been formally constituted as part of the WMO/ICSU World 
Climate Research Programme (WCRP) by the Joint Scientific Committee (JSC) 
for the WCRP and the Committee cn Climatic Changes and the Ocean 
(CCCO), under the guidance of the joint JSGCCCO Working Group on A.r-Sea 
Fluxes (WG-ASF members are listed in Appendix 1). This step was taken in 
view of the requirement in the WCRP, especially the WOCE and TOGA 

,ucu m section 1 . 1 ; Tor continuous global gndded air- 

sea flux fields exploiting the entire range of available input data including, 
in particular, the new sources of ocean satellite data (e g scatterometer 
winds and alt.metr.c information) that are expected to become available 
operationally in the next two or three years 

Although ocean satellites can provide important data on fields at the sea 
surface as input for the flux computations, the derivation of valid geophysi- 
cal data from the sensor signals of these satellites is not straightforward 
Normally, auxiliary meteorological ur oceanic data are needed. The satellite 
data therefore needs to be combined with conventional meteorological and 
sea surface data ir a comprehensive data assimilation system. This combi- 
nation should be performed within the physical constraints of an atmo- 
spheric general circulation model, and the expertise for this activity resides 
at the global weathef centres. 

In addition to ocean satellite data, an important data base for air-sea flux 
computations are conventional in situ ocean surface data, which provide a 
calibration reference for ocean satellites and meteorological analyses The 
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GDAP therefore also includes recommentations for upgrading the 
conventional surface data set obtained from ships, buoys and other sources 

The principal goals of the GDAP may be thus summarized, in accordance 
with the terms of reference of the WG-ASF (cf. Appendix 2), as 

1. the development and implementation of comprehensive data 
assimilation systems for the operational computation of global air-sea 
flux fields at global weather centres and 

2. the upgrading of the in situ air-sea flux data collection system. 

1.1 Data assimilation 


The data required for computing global air-sea flux fields originates 
from a variety of sources (satellites, ships, buoys, weather stations) with 
widely differing characteristics regarding coverage, sampling density, 
accuracy and type of variable measured. The reconstruction from these 
data of optimally estimated flux fields which are dynamically consistent 
with the overall atmospheric circulation requires the application of a 
sophisticated data assimilation system based on a high resolution 
global atmospheric model. In view of the known impact of sea state on 
the microwave sensor signals of oceanographic satellites, and also on 
the air-sea transfer processes themselves, it is furthermore desirable to 
incorporate wave data and a global wave model within the 
assimilation system. 

The data assimilation requirements for climate research ^nd opera- 
tional weather and wave forecasting are essential identical It there- 
fore appears logical to develop a single comprehensive data assimila- 
tion system for both applications. Fig. 1 (from Hasselmann, 1985) shows 
an outline of the proposed structure of such a data assimilation system 
in relation to the ocean satellite receiving stations and off-line data 
processing and archiving facilities. The plan foresees two mam 
extensions of existing operational data assimilation systems the 
incorporation of the sensor algorithms within the data assimilation 
system and the inclusion of a global wave model (GWM). 
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— 9 ‘ — General structure of wind and wave data assimilation system including ocean satellite da 
(from Hasselmann, 1935). 
















A quasi real-time, operational data assimilation system offers a number 
of advantages: 


The system can provide sophisticated continuous data quality 
control, with the possibility of immediate feedback to operational 
centres of satellites or other data collection systems if data errors 
are detected Experience ows that a comprehensive data 
assimilation system encompassing a diverse data set is significantly 
more effective in detecting subtle errors in individual system 
components than quick-look testing of separate sub-systems 

The system can provide continuous statistical evaluation of the 
performance of satellite sensor algorithms, again offering a rapid 
feedback if algorithm improvements are needed. 

The assimilation system can be applied during measurement 
campaigns, e g. in the initial calibration and validation phase of a 
satellite, yielding dete which cen be fed » n nuas* rpal-t»mp 
into the campaign operation, in this way it is possible to 
continually monitor the performance of the satellite or other 
measurement systems during the campaign 

More efficient satellite sensor algorithms can be developed which 
make use of first guess wind and wave fields provided by the 
models This is particularly important if the cross coupling of 
different sensor signals (e g wind and waves) is to be properly 
incorporated in the algorithms. 

Since the analysis system is designed to cope with the complete 
data flux in quar real-time, it must process the data at least at the 
same rate as it is produced. This avoids one of the recurrent 
problems of off-line processing systems: a continuously increasing 
processing delay, or, as the only alternative, the restriction of the 
data analysis to limited time periods 

Although a quasi real-time data assimilation system clearly has 
considerable advantages, a basic shortcoming of an operational system 
is that the homogeneity of the data set is endangered through 
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inevitable changes in the analysis procedure as improved assimilation 
systems and sensor algorithms are developed. For applications in which 
data homogeneity is essential, it may therefore be desirable to 
reanalyze the entire past data set with the latest assimilation system. 
Thus, although the major thrust of the GDAP is directed to an 
operational system, the option of a later re-analysis should be kept 
open. This requires the establishment of efficient archiving and 
retrieval mechanisms for all input data. A re-analysis would clearly be a 
costly undertaking, but this may appear less forbidding when more 
powerful computers become available later in the programme. 

The programme required to develop and implement a comprehensive 
air-sea flux data assimilation system is extensive and can be carried 
through successfully only with wide international collaboration. For 
this reason it is proposed that efforts should be concentrated on the 
implementation of assimilation systems at the two major global 
forecasting centres supported by larger scientific communities, namely 
the European Centre for Medium Range Weather Forecasts (ECMWF) 
and the US National Meteorological Center (NMC). However, efforts by 
individual nations to develop similar systems at their national weather 
centres should also be encouraged, since the intercomparison of 
different products from several centres will provide a broader base for 
assessing the reliability of the methods. 

1 .2 In situ data 

The second task of the Working Group is to recommend and initiate 
measures to improve the data base of in situ reference measurements 
which are needed for the validation of flux estimates inferred from 
satellite data or derived from global meteorological analysis The 
Working Group considers that implementation of improvements in the 
quality, and efficiency of utilization, of existing measurement 
programmes would significantly upgrade the available data base. 
However, where there are potential deficiencies these must be 
identified and, where necessary, the initiation of new scientific field 
programmes to produce in situ reference measurements should be 
encouraged. 


The operational in situ systems, which include moored and drifiting 
buoys, the Voluntary Observing Ships (VOS), and island and land 
meteorological stations, are the responsibility of the World Weather 
Watch (WWW) of WMO. Recognizing that it is of primary importance 
that the in situ observational systems required for initialization of the 
atmospheric GCMs be continued, and where possible, improved, the 
Working Group strongly supports actions aimed at fulfilling the WCRP 
requirements for WWW systems as defined in the First WCRP 
Implementation Plan (WMO, 1985). 

A specific Working Group project, described in Section 4 12, will 
identify data from previous scientific field programmes suitable for use 
in verifying the model derived surface fields. The project will also assess 
the quality of avai ! able in situ data, and determine any additional 
quality control procedures which may be required. The results will 
enable the Working Group to assess the need for further research field 
experiments, and the degree to which, with improvements, the 
existing operational in situ data can be used for verification purposes 
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3. Overview of Work Plan 


The work plan of the GDAP may be divided into two stages: 

a preparatory phase, from 1988 until 1990. and 

a full implementation phase, from 1991 onwards. 

3.1 Preparatory phase. 1988 - 1990 

During this phase a series of investigations need to be carried out on 
methods of assimilating data from different sources, addressing m 
particular the new forms of data which will become available from 
oceanographic satellites. Methods of assimilating single level data (i e 
surface data) will require special attention. Other questions which 
need to be addressed include the dependence of the numerical 
satellite sensor algorithms on the assimilation scheme, the incorpo- 
ration of the algorithms within the assimilation system itself and the 
interdependence of wind and wave data in both the sensor algorithms 
and the assimilation scheme. 

Alternative flux parameterisations need to be tested with regard to 
their impact on oceanographic models (tropical and global ocean 
circulation models and surface wave models). The appropriate 
numerical experiments will be earned out in collaboration with the 
WOCE and TOGA Numerical Experimentation Groups and the Wave 
Modelling Group (WAM). 

In preparation for a possible later reanalysis project (see Section 4 1 3), 
investigations will also need to be carried out on different methods of 
reconstructing improved surface flux fields from past data, including 
techniques which do not necessarily require the full data assimilation 
cycle. 

In addition, work will continue on a number of projects initiated 
before the begin of the preparatory phase of the present programme. 
These include: 
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An ECMWF project requested by the TOGA Scientific Steering 
Group to compute all air-sea flux fields beginning 01 .01.1985, 
applying the present data analysis schemes of ECMWF to the 
Centre's archived data set (see also Section 2.2 and Appendix III). 

A joint WAM-ECMWF project for the development and testing of 
an ocean satellite data assimilation and quality control system 
using SEASAT altimeter, scatterometer and SAR data in 
conjunction with conventional meteorological data 

3.2 Full implementation phase. 1991 onwards 

At the end of the preparatory phase, continuous, operational data 
assimilation systems should be implemented at ECMWF, NMC and 
possibly other centres. The assimilation schemes should provide 
6-hourly flux fields at the resolution of the then operational global 
forecast models (~ 1° x 1°). The assimilation scheme should encompass 
all relevant oceanographic microwave sensor data, including both 
wind and wave information, and should also include the ocean satellite 
sensor algorithms. The full operational phase will be preceded by a test 
phase which could coincide with the calibration/validation phase of 
ERS-1 
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APPENDIX IV 


• Time Table for Global Data Assimilation Programme 
for Air/Sea Fluxes 


Preparatory Phase 


I- 

I Full Implementation 
Phase 


no. project 1988 


1989 


1990 


1991 


1992 


1993 


4.1 


data 

arch i val 


minor mod i - 
fication of 
data sets 
archival 
specifica- 
tions : 
definition 
of wave data 
subsets 


Operational 

implementation 


4.2 Scattero- 
meter 
data 

Development of 
algorithm 
suitable for 
data 


Operational 

implementat'on 


assimi 1 ation . 
incorporation 
of algorithm in 
assimilation 
system 

1 



1 

Test of system 
against SEASAT 
and synthetic 
data 


4.3 Coupled 
wind and 
wave 
sensor 
algo- 
rithms 

Development 
of coupled 
algorithms 
for scattero- 
meter, SAR 
and altimeter 

(Test against 
1 SEASAT data 

Operational 
Imp 1 ementat i on . 
Improvement of 
algorithms 
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4.4 Wave 
data 
assimi- 
lation 


Development of methods 
for assimilating SAR 
altimeter and buoy 
wave data in wave 
model s 

Test of 
techniques 
us i ng 

SEASAT data 

Pre-operational 
test in global 
wave model 


Operational 
imp 1 ementat i on 


4.5 Quality lOevelopment 

control 1 of techniques 
for new 
.ocean 

satellite Test against 

and sur- SEASAT and 

face data existing sur- 

face data 


Operational 

implementation 


Imple- |Pre-o 

menta- | imple 

tion of 

wave dependent 
flux parameteri- 
zations 

(in particular 
surface stress 
and energy flux) 


Pre-operational 

implementation 


A Operational imple- 
mentation of agreed 
upon parameteriza- 
tion 


4.7 Improve Test alternative flux 

sensible parameteri sations 

and (WGNE project and 

latent (operational 

heat and diagnostic) 
fresh (see also project 3.9 

water flux 
parameterisations 
( atmospheric 
model experiments) 


Operational imple- 
mentation of agreed 
upon flux parameter- 
izat ions for 
sensible and latent 
heat and fresh water 
flux. 
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4.8 Assessment 
of opera- 
tional model 
radiation 
flux 

computations 


Comparison 
with results 
of JSC radia- 
tion flux 
group 


A 


4.9 Assess- 
ment 
of air- 
sea 

fluxes 
in ocean 
ci rcul ation 
expe riments 


Numerical model experi- 
ments to contribute to 
(1). (2) and (3) 

(TOGA and WOCE-NEGs) 


4. 10 Assess- 
ment 
of 

stress 
parameteri sation 
in global 

wave mude 1 5 


Quasi-operational 
quality control as 
contribution to (1) 


4. 11 Compu- 
tation 
of SST 


Improvement of SST com- . 
putations using & 
improved satellite sen- | 
sor algorithms and ocean| 
model data assimilation | 
techniques j 


Operational imple- 
mentation of agreed 
upon radiation flux 
parametensation 


Further numerical 
experiments for 
general TOGA and WOCE 
objectives (indirect 
qual ity control ) 


Continual operational 
qual i ty control 


Operational imple- 
mentation of re- 
commendation SST 
analysis method 


4.12 In 
s i tu 
data 

Upgrading of in situ ^ 
data collection system 

Implementation of 
upgraded in situ 
data system 

4.13 Reana- (Assessment of homogeneity 
lysis 
project 

6a 

/ of operational flux 
computat ions 

6b 


/\ : Decision points 

Recommended procedures to be used in operational 
Implementation phase: 

^ Momentum (stress) and (mixing) energy flux 


/2\ Sensible heat, latent and fresh water flux 
Solar and infra-red radiation 
A SST computations 

/$\. in situ data collection and quality control 
system 

i f^eS . ^ Possible begin of flux reanalysis projects 

(start of project and period of past data to be 
analysed to be decided later depending on 
experience with real-time operational system). 
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